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CXI. The Optical Effects of Radiation Induced Atomic Damage in Quartz 


By KE. W. J. Mrroneti and E. G. 8. Paar + 
The University of Reading ¢ 


[Received April 20, 1956] 


SUMMARY 


Measurements of the optical absorption of neutron irradiated crystalline 
and fused quartz have been made in the wavelength range 1450-10 000 & 
(8-5-1-2 ev). Two absorption bands associated with atomic displace- 
ments in crystalline quartz have been found—C at 5-7ev and E at 7-6 ev. 
Optical bleaching of the C band, when illuminated with light of the 
wavelength absorbed, occurs even in the purest specimens. A preferential 
optical and thermal bleaching is observed at 5-4.ev during the early 
stages of the bleaching of neutron irradiated specimens. It is suggested 
that this effect is associated with electron release and capture, without 
intermediate retrapping, by a close complementary pair of traps. The 
model which seems to fit the results best is that the C and E defects are 
oxygen ion vacancies and interstitial oxygen ions, which produce absorp- 
tion bands after trapping electrons and holes respectively. The C and # 
bands are also found in fused quartz. 

Results obtained using thin specimens show that the absorption in 
the long wavelength tail of the absorption edge was enhanced during 
irradiation. The absorption limit in irradiated crystalline quartz is at 
8-1 ev and is the same as that found in both unirradiated and irradiated 
fused quartz specimens. The extension of the tail is interpreted as 
arising from the formation of regions of amorphous material as a result 
of local melting. This is taken into account in estimating the con- 
centration of point defects on the basis of collision theories. All the 
damage is annealed by heating at 950°c. 


§ 1. INTRODUCTION 


THERE is a wide transparent region in quartz between the beginning 
of the infra-red absorption (2-7 microns) and the characteristic electronic 
absorption edge (1450 4). For this reason quartz is a favourable material 
for studying point defects by the measurement of their associated optical 
absorption. 

In spite of being one of the most perfect crystalline materials quartz 
nevertheless contains defects in measurable concentrations. For example, 
See sey so ipeieranits oot ie 

+ Now at Radar Research Establishment, Malvern. 
+ Communicated by the Authors. 
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the smoky colour which is induced in most specimens by ionizing radiation, 
and which occurs naturally in ‘smoky quartz’ has been shown by 
O’Brien and Pryce (1954), Griffiths et al. (1954), Ditchburn et al. (1954),7 
O’Brien (1955) and Mitchell and Paige (1955) to be due to substitutionally 
incorporated aluminium impurity. 

The main interest of the work described in the present paper has been 
to determine those effects which are produced by particle irradiation, 
and which are associated with atomic displacements. This is the only 
practical way of introducing interstitial atoms and vacant sites into 
quartz. To utilize fully the wide transparent region our measurements 
of these effects have extended to the absorption edge. 


The concentration of defects 

Point defects in crystals, either existing from growth or induced 
by irradiation, may give rise to an observable absorption spectrum 
directly. Thus, many impurities in the alkali halides (e.g. Tl, Johnson 
and William 1953) give absorption bands, as also do anion vacancies in 
potassium chloride (« band, Seitz 1954) and vacancies in diamond 
(Clark et al. 1956). These defects do not need to capture or lose electrons 
before becoming observable so that the concentration of absorbing centres 
(WY 4) is equal to the concentration of defects (N,), being related to the 


absorption by PARSON Af ne ae 


where .(cm—?) is the absorption coefficient at the band maximum, 4 E(ev) 
the width of the band at half-absorption and f the oscillator strength 
for the transition involved. 

Alternatively, the point defect may be observable only after capturing 
or losing an electron—for example, the aluminium defect in quartz 
gives rise to absorption after losing an electron. In this case the defects 
may become observable after irradiation with ionizing radiation. The 
saturation value of the concentration of centres will, in general, not be 
equal to the concentration of defects but will depend on the nature of the 
ionizing radiation, the cross-sections for the capture and release of electrons 
by the electron traps, and for the capture and release of electrons by the 
positive hole traps. The concentration of absorbing centres will be 
given by (1) but M,=aNy,, where x is some function of the various 
cross-sections. 

If the particle radiation is producing defects of this kind then it is 
only because all particle irradiation produces ionization that the defects 
lead to observable absorption. Providing that the ionization is sufficient, 
so that at any stage of the irradiation ‘ saturation’ of the electronic 
effects is achieved, the value of “4 may be used to follow the rate of 
formation of defects. This can apply to neutron irradiations since 
ionization processes will occur as a result of the movement of those knocked 
on atoms which are charged, and as a result of the y-irradiation which 
accompanies all pile irradiations. 


{ This paper will be referred to as I, 
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Most absorbing centres in ionic and semi-ionic solids seem to be of 
this type. It is often implicitly assumed, however, that a—1, as in the 
calculations of defect concentrations in fused quartz by Yokota (1952) 
and Levy and Varley (1954), both of whom also assumed that f=1. 


§ 2. MATERIALS AND METHOD 
2.1. The Quartz 

Crystalline quartz is normally transparent in thin sections from 0-46 ev 
(2:7 microns) to 8-3 ev (1500 A) and our measurements indicate that the 
absorption edge, which we take for convenience as the photon energy at 
which the absorption coefficient is 400 cm~1, is at 8-5 ev (14504). At 
this energy the absorption coefficient is increasing at a rate of about 
200 cm~? per 0-lev. That the absorption edge defined in this way is 
a property of pure crystalline quartz is suggested by the fact that the 
same value was obtained for both synthetic and natural crystals. It is 
unlikely that this would have been the case had we been working in the 
tail of the absorption edge, whose form can be influenced by crystalline 
imperfections. In the fused quartz used in the present work the 
transparent region extended to 8-1 ev (14504). The absorption below 
8-5 ev in fused quartz can be regarded as an enhanced tail of the crystalline 
quartz absorption edge arising because of the absence of long range order. 

The crystalline quartz used in the present work was either synthetic, 
grown by the G.E.C. Research Laboratories, Wembley, or naturally 
occurring Braxilian quartz. The fused quartz was of optical grade 1 
obtained from the Thermal Syndicate. 


2.2. The Optical Measurements 

In the wavelength range 2100 to 10 000 A transmission measurements 
were made using a Hilger UVISPEK spectrophotometer. From 1450 
to 22004 the transmission was determined photographically using 
a fluorite vacuum spectrograph of the type described by Cario and 
Schmitt-Ott (1931). The arrangement of the apparatus is shown in 
fig. 1. Exposures were made on Ilford Q1 plates, the wavelength calibra- 
tion being obtained from the hydrogen emission lines from the source, 
and the absorption lines of oxygen which was deliberately introduced 
for this purpose during one exposure. The hydrogen lamp was of the 
type described by Heddle (1952). The grids, which were inserted into 
the beam to determine the characteristic curve of each plate, were similar 
to those used by Ditchburn and Heddle (1953) but in the present case 
they could be calibrated in situ, using a crystal of irradiated quartz whose 
transmission between 2100 and 2200 4 had previously been determined 
in the spectrophotometer. Specimens were mounted in the holder 
illustrated in fig. 2. The crystal is held in the vice V which could be 
moved by the screw driver D, controlled from the outside through a 
Wilson seal, so that the specimen could be in, or out of,the beam. Low tem- 
perature measurements were made with the well (W) filled with liquid air, 
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Fig. 1 
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The arrangement of the apparatus for the vacuum ultra-violet measurements : 
L, water cooled, d.c. hydrogen discharge lamp; E, exposure device ; 
S, specimen holder ; G, grids for the relative intensity calibration of the 


plates. 
Fig. 2 
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Diagram of specimen holder for use with the vacuum spectrograph in the 
liquid air to room temperature range. 
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By these two methods the percentage of incident light transmitted (7’) 
by a thickness d cm of crystal has been determined from 1450 to 10 000 4. 
The measurements were made under conditions of normal incidence so 
that the absorption coefficient (1 em~1) is defined by 


DET exp (apd) 7 eee G8) 


where 7’, is the percentage of incident light transmitted by a non-absorbing 
specimen having the same refractive index. If there is no scattering of 
light from the two surfaces and if, as in the present case, multiple 
reflections can be neglected, 7') can be expressed in terms of the refractive 
index (n) by me a\2 
T,=(1—R)? “where R= ) sees eas ( 3) 


Fig. 3 


WAVELENGTH (A) 


2000 _ 1800 1600 2000. 1800 1600 
i ae I 1s eel ia Pe ae 


MS CRYSTALLINE FUSED 


| | = 
3 7 8 


6 7 8 
PHOTON ENERGY (eV) 


Values of 7’, (eqn. (2)) as a function of wavelength in the vacuum ultra-violet 
region for crystalline and fused quartz. 


Between 2000 and 100004 the values of n for crystalline and fused 
quartz were taken from the International Critical Tables (1933). In 
crystalline quartz the difference between 7’ calculated from n for the 
ordinary ray and n for the extraordinary ray was less than 0-2%. Such 
a difference is not significant in the present work so that values of 7) 
calculated from n for the ordinary ray have been used. In the region 
2000-1450 A the reflection coefficients determined by Johnson (1941) 
have been used for fused quartz. For crystalline quartz we have deter- 
mined 7',(=82%) at 16104 from transmission measurements on two 
unirradiated crystals of thickness 0-01 and 0-2 cm. We have extrapolated 
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the data in the tables, which goes down to 1850 4, through the 1610 4 
point, giving the curve of 7’, as a function of wavelength as shown in fig. 3. 

In using values of 7', determined from unirradiated specimens to 
calculate the induced absorption coefficients it is assumed that: (a) the 
irradiation does not significantly affect n, (b) the surface topography is 
not changed sufficiently to cause significant scattering of light. 

Some change of n must be expected to be associated with any absorbing 
centre, but for absorption coefficients up to the highest measured in the 
present work (500 cm-1) the consequent change in reflectivity can be 
neglected. In order to check the second assumption the transmission 
in the visible was measured for a synthetic specimen which had received 
a large neutron dose (1019 n° cm-*). The weak Al type band present in 
the crystal was completely radiation bleached and the measured trans- 
mission was within 1%, of the value of 7’, (eqn. (2)) calculated from the 
index data. Thus no appreciable surface scattering was introduced 
by this dose. 


2.3. The Irradiation 


The specimens have been irradiated in the Harwell Pile (BEPO) in 
the centre of a cell of the graphite lattice at an ambient temperature of 
about 50°c. The neutron spectrum in this position is given by 


n(e)de=A (l/eg?. « . exp (—e/eg)+0-05/e)de. . . . « (4) 


where €, is about 0-025 ev, n(«)de is the flux in neutrons per cm? per sec in 
the energy range « to e+-de, and A is a constant depending on the power 
level of the pile and is given approximately for each irradiation. 

The second term does not apply in the range of energy up to 0-1 ev 
so that we can write the integrated thermal flux (D,;) as 


O-1 
D,= | ne) de=0-91 A. 
0 


The flux of high energy neutrons is then given by 
n(e)dersO-055D,defe . . . . . . . . (5) 


provided that «>100 ev. 

Neutrons will not be expected to produce displacements of silicon 
or oxygen atoms when the neutron energy is less than about 100 ev so 
that (5) may be used to determine the energies of the primary knock 
on atoms. The neutron doses quoted in this paper are values of the 
product (D, x exposure time in seconds)=D. 

The induced radio activity was allowed to decay. In some cases 
measurements were made and they indicated that the activity was 
associated with traces of polishing powder which had adhered to the 
ground faces perpendicular to the polished surfaces. Autoradiographs 
confirmed this distribution. | 
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§ 3. EXPERIMENTAL RESULTS 
3.1. The C Band 


In good synthetic or natural crystals in which the concentration of 
aluminium impurity is of the order of, or less than, 0:01 atomic % the 
form of the damage induced CO band (5-7 ev) can be seen. This is shown by 
curve I of fig. 4, which is the absorption curve of a synthetic crystal 
after a dose (D) of 7-4 10!’° cm-?. It will be noticed that even for 
this crystal there is some absorption in the region between 1 and 3 ev 
and the amount is consistent with the aluminium concentration found by 
Brown and Thomas (private communication) and Griffiths et al. (1954) 
in similar synthetic samples and also by Volger et al. (1955) in synthetic 
quartz from Bell Telephone Laboratories. Figure 4 also shows curves 
(II and III) for less pure specimens which had received the same dose. 


Fig. 4 
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Absorption spectra of neutron irradiated crystalline quartz : dose, D, 
7-4 1017 n° cm; curve I synthetic ; curves II and III natural Brazilian. 


These curves indicate that some absorption in the C' region, in addition 
to that induced by the damage, is present and that it is greater for 
specimens having greater A,(2-0 ev) and A,(2:7,ev) band absorption— 
ie. greater aluminium concentration. We have found this to be true 
for all specimens. 

Absorption in the C region can be induced by x-rays. It was shown 
in I that this absorption saturated, so that some of the defects which 
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after ionization give rise to the C absorption must be present in the 
crystals initially. The level of C absorption due to trapping in existing 
defects seldom exceeds 20 cm~! and is usually very much less than this. © 
In a good synthetic or natural crystal it is about 5cm~1. Thus for the 
higher neutron doses the C band absorption is mainly due to displacements 
induced by the irradiation. 

Figure 5 shows the neutron rate curve for the C band in different 
crystals for doses up to 10!9n°cm-*. It is interesting to note that 
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Rate curves for the C band induced by neutron irradiation. 


the points for all specimens eventually lie on the same curve which we. 
call the Damage Rate Curve. For example, the results for the specimen 
whose rate curve is denoted by crosses lie on the Damage Rate Curve 
after a dose of 30x10! n° cm-*. This specimen, judging by the A 
absorption coefficient, was the least pure and it appears that the initial 
contribution to the C absorption has decayed with prolonged dose in 
a qualitatively similar manner to the 4, A, bands, fig. 6. 
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The Damage Rate Curve is slightly curved but the absorption 
coofiicient at the band maximum increases approximately as 10-17 x D. 
The half-width of the band is close to 1-4 ev so that 


(oped eee 4 10st XD 10-18 ef, 
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Here /, is the concentration of C centres and f, the oscillator strength 
for the transition. 
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Rate curves for the A, band induced by neutron irradiation. 


3.2. Absorption at Photon Energies Greater than 5-8 eV 


The extent to which defects initially present in the crystals give rise 
to absorption in the vacuum ultraviolet after irradiation with ionizing 
radiation is shown in fig. 7. There is an additional band (D band) 
having a maximum absorption at 7-1; ev of 7 em~! for this crystal after 
a 125 Mr dose of 50 kv x-rays. The values of absorption coefficient for 
energies greater than 7:4 ev are uncertain because of the subtraction of 
the absorption in the tail of the absorption edge which, owing to the 
thickness, was significant in this specimen. 

The effects of neutron irradiation are illustrated in figs. 8 and 9. 
Curve 1 in fig. 8 refers to a synthetic specimen which had received a dose 


Fig. 7 
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The absorption spectra of x-irradiated natural crystalline quartz: I, - — — — - 
UVISPEK measurements, —————— vacuum ultra-violet measurements. 


The absorption in the unirradiated crystal, curve II, has been subtracted 
from the measured absorption to give curve I. 


Fig. 8 
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Vacuum ultra-violet measurements of neutron irradiated crystalline quartz : 
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of 75x 10" n° cm-*. It will be seen that the C band peak is also found 
in the vacuum ultraviolet measurements. However, the transmission 
at 6-9 ev had fallen to less than 3% so that values of absorption coefficients 
were not obtained for higher photon energies with this specimen. 

For a thinner specimen (0-011 cm) which had received a neutron dose 
of 17X10!" n° cm-?, the spectrum shown in fig. 9 was obtained. The 
most prominent feature is the H band which has an absorption peak 
at about 7-6 ev and a half-width of 0-5 ev. After the decrease in absorp- 
tion on the high energy side of 7-6 ev the absorption rises rapidly to 


Fig. 9 
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A complete absorption curve for neutron irradiated crystalline quartz. 


immeasurably high values, the last measurable value being 300 cm™ 
at 8-lev. No evidence has been found for a decrease in absorption 
between 8-1 and 8-5ev, the induced absorption limit being similar to 
that found in both unirradiated and irradiated fused quartz at 8-1 ev. 
If the possible absorption tail denoted by the broken line in fig. 9 is 
subtracted from the measured curve the H band appears as a symmetrical 
band. At —195°c the # band is sharper and the peak has moved to 
a higher energy (fig. 8). 
Because of the low transmission at the H band peak the absorption 
coefficient is uncertain within +20%. However, within these limits 
the absorption due to the EZ band has been found to vary linearly with 
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It has not been possible to discover whether the D band is enhanced 
by neutron irradiation because of the overlapping of the strong C and # 
bands in this region. If it is enhanced it is clear that the rate of increase 
is very much less than that of either the C or # bands. 


3.3. Optical Bleaching of the Induced Absorption 


Induced absorption bands which decrease in strength during 
illumination must be associated with defects which have trapped electrons 
or holes. The work of Choong (1945) and that reported in I has shown 
that in quartz both the A and C bands decrease upon illumination with 
light in the C band. The bleaching rates are, however, small so that 
it is not feasible to carry out the experiments using monochromatic 
light. 

“Tllaminstion with light in the A bands leads only to a small optical 
bleaching effect (table 1) so that for experiments in which bleaching 


Table 1. Effect of Iumination for 60 hr with Light of Photon Energy 
less than 3 ev 


(High pressure Hg lamp and filters) 


Transmission Before After Instrumental 
measured at Illumination Illumination 


2-7 ev 
5-7 ev 


rates have been measured using light between 2100 A and 10 000 A the 
main effect is due to absorption in the C region. 

In a y-ray (Co) irradiated specimen in which the dose was such 
(60 Mr) that only ionization effects were being studied the effect of optical 
bleaching is shown in fig. 10. The difference curve indicates that the 
spectrum bleaches uniformly. This is in contrast to the optical bleaching 
of neutron irradiated specimens reported in I. In this case, under the | 
same conditions of illumination, preferential bleaching occurred in the 
region of 5-4 ev (fig. 15 of I). The preferential bleaching has now been 
shown to occur during the initial rapid bleaching of the C band shown in 
fig. 16 of I. This may be seen more clearly in fig. 11 where three difference 
curves are plotted for different parts of an optical bleaching rate curve. 
After 80 minutes bleaching there is a pronounced maximum at 5:4 ev 
which becomes less pronounced during the next 180 minutes, and is 
absent from the difference curve obtained after a further 320 minutes. 

Experiments which have been carried out on a neutron irradiated 
(92x 107 n° cm™) synthetic crystal indicate that the presence of 


Fig. 10 
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The absorption spectra of y irradiated crystalline quartz: curve I after the 
irradiation ; curve II after subsequent ultra-violet illumination for 22 hr ; 
curve III the difference spectrum. | 
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The optical bleaching of a neutron irradiated crystal by illumination with light 
in the A and C bands: curve I, the change after the first 80 min ; 
curve II, after the next 160 min ; curve III, after the next 320 min. 
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the A bands is not essential for optical bleaching. These curves (fig. 12) 
show that the C band is bleached and that during the first 30 hours there 
is preferential bleaching in the region of 5-4ev. (In the specimen which 
had received the lower dose 2:4 1017 n° cm above, the preferential 
bleaching effect was absent after 5 hours). During the following 80 hours 
bleaching is less rapid and the preferential bleaching is relatively much 
smaller. Presumably one of the vacuum ultra-violet bands is also being 
bleached but optical bleaching measurements on these bands have not 
yet been made. 

If the illumination is with light of all wavelengths up to the absorption 
edge it is found that the C band is bleached much more rapidly in relation 


Fig. 12 
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The optical bleaching of neutron irradiated synthetic quartz: dose 
92 «1017 n° em-?; curve I, no illumination ; curve II, 30 hr; curve III, 
after a further 80 hr. 


to the A band than if illuminated with light in the 2100-10 000 A region. 
This experiment was carried out using a neutron irradiated natural 
crystal and is another indication that some of the vacuum ultra-violet 
absorption must be optically bleached. 


3.4. The Effects of Heat Treatment 


It is necessary to distinguish two processes which can lead to decreases 
in absorption—the thermal activation of electrons, and the thermal 
activation of atoms or vacancies. For convenience we refer to the 
former as ‘ bleaching ’ and the latter as ‘ annealing ’. 

The A and C bands may be bleached at temperatures less than 500°c 
and it will be shown later that at this temperature there is little annealing 
taking place. The isothermal bleaching curves do not follow a simple 
rate law but there are two interesting features. Firstly, in neutron 
irradiated specimens there is preferential thermal bleaching in the region 
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of 5-4 ev which is not found for y- or X-ray irradiated specimens. This 
is shown by the curves of fig. 13. Secondly, it is possible to deduce 
activation energies for the electron excitation from successive isothermal 


Fig. 13 
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hing of radiation induced absorption in crystalline quartz. 
oe erates ates the change in absorption after heating for 9 hr at 
240°c for: I, x-irradiated crystal; I, neutron irradiated crystal 
(small dose), = —-—- = after the first 10 min and- — -— -—- - after 

the next 8 hr 50 min; III, neutron irradiated crystal (large dose). 


bleaching curves on the same specimen. This procedure, e.g. Overhauser 
(1953), is independent of the rate kinetics. The energy ¢<’ is given by 


(2), /@d-mele-2) 0 
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For crystalline quartz, measurements at 215 and 335°c on one specimen 
and also at 220 and 260°c on another both gave «’/(2) equal to 0-75 ev. 


Annealing Experiments 

High temperature heat treatments were carried out in either air 
or vacuum, the crystal being contained in either a platinum crucible 
or wrapped in platinum gauze. From 950°c it was necessary to cool 
the crystal slowly to prevent cracking and in these experiments specimens 
were cooled at a controlled rate (approx. 40°c/hour) for 24 hours. To 
decide the extent to which the atomic defects are annealed at any given 
temperature it is necessary to re-irradiate the specimen after the heat 
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The effect of a neutron irradiation-heat treatment-neutron irradiation cycle 
on the absorption spectra of natural crystalline quartz: curve I, after 
the first irradiation ; curve II, after the second irradiation ; curve III, 
after the same irradiations but with heat treatment at 950°c for 24 hr 
in air between irradiations. 


treatment. Ideally x-ray irradiations to saturation are desirable for 
this purpose but in many cases this process is so long as to make the | 
experiment impossibly slow. The alternative of a standard x-ray dose 
for the second irradiation cannot be relied upon since Forman (1951) 
has shown that heating at even 400°c can affect the rate of visible colour- 
ing although the saturation value appears unaltered. In the present 
experiments we have used a second pile irradiation to determine the 
effect of the heat treatment, this being done in conjunction with an 
irradiated, but unheated, control specimen, 
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Three specimens were cut from the same crystal, all with the c-axis 
perpendicular to the polished surfaces. These specimens were placed 
in the same container and received a dose of 241017 n° cem-2. After 
measuring the absorption spectra between 1-2 and 5-8 ev, one specimen 
(N32) was heated for 45 minutes at 500°c in air, another (N34) was 
heated in air for 24 hours at 950°c and the third specimen (N33) was 


untreated. All specimens were re-irradiated (the C absorption of the 


control increasing by 8-8 cm!) and the absorption spectra redetermined. 
Curve I of fig. 14 is the absorption spectrum of one of the crystals after 
the first irradiation. This was typical of all three specimens, the 
differences being +0-05 cm-! for uw, and +2 cm! for u,. Curve II is 
the absorption spectrum of N32 after the second irradiation and the 
absorption spectrum of N33 was similar. The major change shown 
by the second irradiation was the increase in the absorption of the A bands 
of N34 after the high temperature heat treatment (curve III). The same 
large increase in A band absorption was found for N33 following an 
irradiation after high temperature heat treatment in vacuum. Thus ~ 
the effect is not significantly sensitive to the ambient atmosphere. Finally, 
heating N34 for a further 24 hours in air after the second irradiation 
produced a further enhancement of the A bands, which was less, however, 
than the original increase. 

In view of the evidence that the A bands are associated with 
substitutional Al, the first conclusion to be drawn from these experiments 
is that heating to 950°c after neutron irradiation leads to a higher con- 
centration of substitutional Al. This in turn indicates that normally 
some of the Al impurity is present in interstitial positions.— Then 
during heat treatment the interstitial Al atoms diffuse at a similar rate 
to the interstitial Si atoms produced by displacements, so that some of 
the Si vacancies are filled by Al atoms. 

It has already been pointed out that the C band absorption due to 
defects initially present in the crystal seems to be related to the A band 
absorption, that is to the substituted Al concentration. The heat 
treatment results lead to a similar conclusion, for an increase in C' absorp- 
tion always accompanied the A band enhancement in spite of whatever 
annealing processes were taking place (cf. curve IIT, fig. 14). 

To examine this quantitatively the heat treatment results have been 
analysed using the two alternative relations : 


(a) o=e, -hD ake en aul Eine es) 
(b) po=2’'tatBD+Dy be Pope iLO) 


where «, «’ and f are constants, D is the neutron dose and. Dy is a 
contribution to the absorption from defects in the crystal prior to 
irradiation, for example, the tail of the D band. It is clear that using (9) 
« and B cannot both be constants applying to all three specimens for 


+ The non-linear increase in ha with Al concentration reported by Ditchburn 
et al. (1954) is a further indication of the presence of interstitial Al 
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after the first irradiation zg was not constant. We assume that f, the 
damage component per unit dose, was constant. In the comparison of 
results in table 2 we consider the extremes of no annealing and complete 
annealing for each temperature. 

There is not a great deal to choose numerically between the two 
analyses, what difference there is being in favour of (10) in which «’ is 
constant for all quartz crystals. In either case the results lead to the 
conclusions that there is very little annealing at 500°c and that annealing 
is complete at 950°c. These conclusions are also supported by the 


Table 2. Analysis of Heat Treatment Results 


: Heat FTP 9s en : p,.cm~* | , Clie 
Sy ee treatment | observed ee calc. 9 | calc. 10 
N 32 No annealing 30 30 
Mi aie 500°c in air | 30+1 
tee Complete a : 
Dore annealing a a 
No annealing 35 33 
950°c in air | 30-5+1 — 
N34 Complete 30 28 
gael by annealing 
a’ = 1 % 3 5 =. 
5 Be ANY: No annealing 47 46 
950°c in air | 32+1 tit 
Complete 34 99 
annealing 
N33 No annealing 51 48 
oil 950° in : 2s _ 3 
a eeisd. vacuum saan) Complete 38 34 
Dg=2-6 annealing 


observations on the absorption edge. The same result has been found 
by Wittels (1953) for the annealing of the density changes, while Berman 
(1951) has found that the irradiation induced thermal resistance was not 
completely annealed at 700°c. 

Measurements in the vacuum ultra-violet were also made on the— 
specimens N 33/34. It was found that after heating at 500°c for 2 hr 
the # band was still present, although reduced in strength (from 100 
to 35cem~'). After heating at 950°c the H band was completely bleached 
and the absorption in the region of the edge was identical with an 
unirradiated crystal. 

While the above heating results refer to high temperature heat 
treatments carried out after neutron irradiations, it has been found that 
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the amount of A band absorption induced may also be altered by heat 
treatment before irradiation. Two specimens cut from the same block 
were irradiated, one of them having been previously heated in air at 950°C 
for 24 hr while the other was untreated. The value of the absorption 
coefficient at the A, band peak was 17° higher in the heated specimen. 
This is much larger than the variations in the A, band absorption amongst 
non-banded specimens cut from the same block (--2%), although much 
smaller than when the specimen has undergone an irradiation-heat- 
irradiation cycle. We conclude that it is possible by heat treatment to 
increase the concentration of substituted Al in an unirradiated crystal. 
It must be assumed that there is a variable concentration of Si vacancies 
existing from growth in unirradiated crystals, for it is improbable that 
sufficient vacancies can be created by heat treatment at 950°c. 


3.5. Observations on Fused Quartz 


As reported previously in paper I for the spectral range 2100 to 
10 000 A, the behaviour of fused quartz is similar to that of crystalline 
quartz and this similarity of the shape, positions and properties of the 
bands has now been found to extend into the vacuum ultra-violet region. 
The local environment of such centres must be similar in the two solids, 
a view which is consistent with the model of fused quartz as a solid 
lacking the long range order of the crystalline form. Recently Weeks 
(1956) has shown that similar magnetic resonance lines are associated 
with the radiation induced damage in the two solids. The optical effects 
are summarized below and some of them are illustrated in fig. 15. 


(a) Corresponding to the A, and A, bands in crystalline quartz there 
isa band (A) with a peak at 2:3ev. This absorption also decreases with 
prolonged neutron irradiation. 

(b) There is a band at 4:1 ev (B,) which is strong after x-irradiation, 
although in some specimens the absorption passes through a maximum 
before reaching a steady value. With neutron irradiation it is bleached, 
in all specimens, more rapidly than the A band. 

(c) There is a band induced at 5-1 ev (B,) in neutron irradiated, but 
not in x-irradiated specimens. The same result was found by Levy 
(1955) but the band was not observed by Levy and Varley (1954). 

(d) A band (C) occurs at 5-8ev and increases with neutron dose, 
but more rapidly than the C band in crystalline quartz. 

(c) In x-ray irradiated fused quartz a band is found at 7:2ev (D) 
as in x-ray irradiated crystalline quartz. It is similarly not known 
whether the band is present in neutron irradiated specimens because of 
the rapid growth of the C and # bands. 

(f) The # band occurs at 7-6ev. In fused quartz, however, the 
E band has been found in a specimen which was irradiated with x-rays 
with a dose 7X that required for the A band saturation. It was not 
observed at saturation. This suggests that in unirradiated fused quartz 
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E defects either exist or can be produced with a small probability by 
x-ray photo electrons. The latter alternative may arise since in fused 


quartz, because of the irregular packing of the SiO, groups, there will © 


exist atoms having binding energies less than in crystalline quartz. 
The band has also been found in y(1 x 10°r) irradiated specimens. 

(y) The absorption limit is found at about 8-1 ev in unirradiated and 
irradiated fused quartz, and is not affected by heat treatment. 


Fig. 15 
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The radiation induced absorption spectra in fused quartz: curve I, after — 
neutron irradiation; curve II, the same specimen after heating at 
520°c showing the presence of the # band; curve III, an x-ray dose ~ 


corresponding to saturation of the A band ; curve IV, the same specimen 
after <7 the saturation dose. 


Since the visible absorption does not increase with neutron dose 
the A band is probably associated with impurity. It has been reported 
by Levy (1955) and also Cohen (1955), that specially purified fused silica 
manufactured by the Corning Co. does not colour with x-irradiation, 
and shows only a small visible colour after neutron irradiation. The 
effect may not be due simply to purification, however, as Yokota (1952) 
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has shown that in fused quartz prepared under strongly oxidizing con- 
ditions the A band is absent. Similarly we have found that the strength 
of the A band induced by a small neutron bombardment, while practically 
unaffected by prior heat treatment at 950°c for 24 hr in vacuum, is less 
if the treatment is in air (table 3). Levy and Varley (1954) have suggested 
that boron is the relevant impurity and that it exists in amorphous 
quartz as sodium tetraborate, the A and C centres being induced by 
hole and electron trapping in the borate. The general similarity with 
crystalline quartz suggests that boron may play a role in fused quartz 


Table 3. The Effect of Previous Heat Treatment on the Absorption 
Induced in Fused Quartz by a Neutron Dose of 4:2 x 10!7 n° em 


Absorption coefficient (em~') 


Specimen Treatment : 
A B, B, C 
F7 os 8:5 71 16:8 20-0 
Fs 950°o 24 hr in vacuum 8:2 7-0 17:3 20-0 
F9 950°o 24 hr in air 5:7 6-5 ies 19-4. 


similar to aluminium in crystalline specimens and exist in place of silicon 
in the SiO, network rather than occur as sodium tetraborate. In six 
of our specimens which were analysed, however, the boron concentration 
was less than 0-0001°% by weight, too small to give an observable 
absorption band by the same mechanism as aluminium, while the 
aluminium concentrations were all around 0-05%. 

As in crystalline quartz, all the defects introduced by neutron irradiation 
may be annealed by heating at 950°c. Thermal bleaching of the A and C' 
bands takes place with an activation energy ¢'/(2) of 1-6ev. A detailed 
account of the results on fused quartz will be published separately. 


§ 4, Discussion oF RESULTS 


The three primary optical effects of the neutron irradiation of quartz 
with doses up to 1019 n° cm~* are: (i) the growth of the C band, (ii) the 
formation and growth of the # band and (iii) the edge shift. The first 
two of these occur at practically the same energies in both fused and 
crystalline quartz, while the edge shift is only found in crystalline 
specimens. All three types of defect are completely annealed by heating 
at 950°c but are hardly affected at 500°c. Neutron irradiation also 
decreases the strength of the visible absorption (A,, A,) bands associated 
with Al impurity—the effect which we have referred to as radiation 
bleaching. 

4.1. Radiation Bleaching 

Continuous radiation bleaching of the A bands occurs only with neutron 

irradiation and we examine first the rates to be expected for the destruction 
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of the defect. There seems to be two possibilities: (a) transmutation of 
Al by thermal neutrons, and (b) displacement of the Al or some of its 
neighbouring atoms. Both these processes, however, give rates which — 
are too low to account for the observed changes. In the case of (a) 
this may be seen immediately from the low cross section for the trans- 
mutation of Al (0-008 barns) while the calculations in (iii) below show 
that a given Si (or Al) atom has a chance of ca. 1/1000 of being displaced 
after a dose of 1018 2° em~-?. Experimentally the decrease is by about 
1/5 after this dose. Even if it is supposed that the displacement of any 
one of 10 surrounding atoms destroys the defect the rate of destruction 


is still too small. 


Fig. 16 
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The absorption spectrum of neutron irradiated crystalline quartz, curve I; 
after optical bleaching, curve II ; after optical bleaching and subsequent 
X-irradiation, curve IIT. 


It is more probable that the effect is electronic. There is no evidence 
of any new centre being formed at the same time as the A bands decrease, 
The following experiment, however, indicates that the ratio (W,/N,) 
for a given type of radiation depends on the relative occupancy of the 
other trapping centres. Such a dependence could provide a basis for 
a mechanism of radiation bleaching. 

A specimen which had been partially radiation bleached (3-0 x 1018 n° 
cm~*) was irradiated with x-rays and no increase of the A,,A, band absorp- 
tion was found. The absorption curve of a similarly radiation bleached 
specimen is shown in fig. 16 (1). This specimen was further bleached by 
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illumination with C (and A) band light, giving the absorption curve IT. 
An x-irradiation of the same dose as for the other specimen led to the 
absorption curve IIT. Since the action of the light cannot be to anneal 
atomic defects the increase of A,, A, band absorption above the original 
level must be due to change in the ratio W,/N,. 


4.2. The Edge Shift 


On the high energy side of the # band the absorption rises rapidly 
at 8-1 ev, compared with the corresponding rise at 8-5 ev in unirradiated 
quartz. The increase at 8-lev is similar to the increase found in 
unirradiated fused quartz. 

In fused quartz the SiO, groups are not regularly arranged. The 
variations in Si—O distances presumably lead to allowed energy levels 
occurring below the energy at which the conduction band is found in 
crystalline quartz, and above the valence band. The absorption spectra 
indicate that the neutron irradiation of crystalline quartz produces levels 
in the same energy range. We infer that some regions of the crystal 
have been converted to a phase analogous to fused quartz. 

It is not possible to estimate the total volume occupied (after a given 
irradiation) by these regions from the present measurements since this 
involves knowing the whole absorption curve up to values of 10°-108 cem—} 
for both the irradiated and unirradiated specimens. Single crystal 
specimens have to be used and for reliable measurements thicknesses 
of 10-5 and 10-*cm are required. As a method of detection of a 
disordered structure, however, the sensitivity is large. The present results 
could correspond to the conversion of only 10% of the total volume of 
the crystal, equivalent to a density change of 0-015%%. 

The occurrence of disordered regions (M regions) for the doses used in 
our measurements (up to 101° n° cm-*) is of some interest following the 
work of Wittels and Sherril (1954), and the discussion of this work by 
Kinchin and Pease (1955). Wittels and Sherril found that after doses 
of about 102° n° cm-2, an amorphous, optically isotropic material was 
obtained having a density 15% less than «-quartz and 3% greater than 
a typical value for fused quartz. For doses up to 7x 10"° n° cm™ the 
damage, as indicated by density and x-ray measurements, could be 
annealed by heating at 930°c, but beyond this dose the heat treatment 
gave polycrystalline «-quartz. The simplest explanation of the trans- 
formation to the amorphous state is that near the ends of the tracks of 
the knocked on atoms collisions are so frequent that we may consider 
a small region being heated momentarily to a high temperature. This 
is the ‘displacement spike’ discussed by Brinkman (1954). In the 
case of quartz it is suggested that the rapid heating and cooling leads 
to the formation of a volume v of amorphous solid. This process would 
lead to a saturation according to a rate 


Vo Vil—exp(—Dav)) se gt ie ag 
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where V,, is the total volume of amorphous solid, V is the volume of the 
crystal, v the average volume of an amorphous region and « is the number 


of knocked on atoms, per unit dose, which produce such regions. 


However, the variation of density with neutron dose observed by Wittels 
and Sherril was not of this form. The relation (11) predicts that the 
slope should decrease with dose over the whole range, contrary to what 
was observed. 

Kinchin and Pease (1955) have suggested that the quartz structure 
becomes unstable when the crystal contains a sufficient concentration of 
point defects and that this is the reason for the rapid increase observed 
by Wittels and Sherril for doses ca. 10?°n° em-*. The smaller decrease 
in density at lower doses is then assumed to be due to the expansion 
produced by interstitial atoms, and not to a smaller number of regions of 
amorphous solid. This is in contradiction to the observed edge shift. 

The experimental results require that regions of amorphous quartz 
are formed at all doses. At the boundaries of such regions there must 
be large stresses, arising because the amorphous material normally 


/ 


occupies a volume about 15% greater than the crystalline material © 


from which it is formed. Fry (private communication) has suggested 
that as the concentration of amorphous regions increases an intermediate 
volume can become sufficiently stressed for the next displacement spike 
to lead to a larger volume of amorphous solid than would be produced 
at lower doses, the additional thermal energy arising because of the release 
of strain energy. 

Suppose v=v,+cD", where c and n are constants, then the rate law 
becomes 


V2=V{1— exp] —aD wt =25 -D" \ oO 3 


To obtain a curve resembling that obtained by Wittels and Sherril we 


have to assume further than when V,—V expansion can take place by — 


the continued production of point defects. Equation (12) can be made to 
fit the three pomts given by Wittels in the range (5-6-6) x 1019 n° em-?, 
but there are no points in the most sensitive range immediately above 
6:6 X 107° 2° om-*. 


4.3. Collision Theory Estimates of the Concentration of Point Defects 


The M regions are not considered to contain significant concentrations 
of point defects, as discussed by Brinkman (1954). Calculations of point 
defect concentrations, therefore, have to allow for the loss of energy of 
the primary knock-on atoms by this process. 

In the case of quartz the primary knock-ons have energies up to 
2-6 x 10° ev for silicon and 4-4 10° ev for oxygen corresponding to the 
highest energy (2 Mev) of the pile neutrons. We have to consider the 
energy loss by four processes: (a) excitation of the electron system 
of the solid, (6) elastic collisions in which the recoil energy is insufficient 
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to produce displacements, (c) elastic collisions which produce displace- 
ments and, (d) the formation of M regions. The loss of energy of primary 
knock-ons by processes (6), (c) and (d) has been discussed by Brinkman 
(1954) and by processes (a), (b) and (c) by Kinchin and Pease (1955). 
We have to apply their results to the spectrum of primary knock-ons 
produced during the pile irradiation. 

The method used is to calculate the upper limit of the concentration 
of defects assuming that no energy is lost by processes (a) and (d). We 
then estimate the maximum effects of (a) and (d) which enables an 
estimate to be made of the lower limit to the concentration of defects. 

The maximum concentration can readily be determined. Neutrons 
in the energy range « to «+de, n(e)de, produce primaries having energies 
E in the range 0 to Fy,ye, where Fy, is the largest fraction of neutron 
energy which can be transferred in an elastic collision. Since the cross 
sections for neutron collisions with oxygen and silicon atoms are isotropic 
and independent of energy in the present energy range, the energy 
density of primaries is independent of energy up to Lyy(—F yy - «yy) 


and is given by (N 0a) m(€)de 


N,dk ae 1 Sete Bee ne (13) 
where NV, is the concentration of atoms a and o, the cross section for n° 
collisions with a. Following the discussion by Kinchin and Pease (1955) 
a primary produces £/2#, displacements in slowing down to rest, as 
long as H<L, (in their nomenclature). This condition is satisfied for 
all the knock-on energies occurring in the present work. Thus the 
maximum concentration of defects will be given by the sum of two terms, 
representing the effects of oxygen and silicon primaries, of the form 


eMx Fyuxt H emx Mica 
pees | | op Nn db+| N, a 
2Eg/F ux 2Eq ao S Eq/Fux Eg 
which, using eqn. (1) and since Pyyeyx >2Hy gives 
Nyx=N aa) Se eer (14) 


We suppose that the knocked on atoms will cease to lose energy by 
ionization when their energy has fallen to H;, the value of which is 
discussed below. The maximum reduction in Nyy will occur when it is 
assumed that above £, all the energy is lost by ionization. Thus taking 
ionization losses into account, the lower limit to the concentration of 


defects is 
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When the energy of knocked on atoms falls to Hy,, the upper limit for 
the formation of M regions as discussed below, we assume that no 
permanent point defects are produced. Thus, considering processes 
(b), (c) and (d) we obtain a lower limit 


€ F € € E 
aa ph rye oN, di {ee ar 
Ey /F ux Ey 24 2 d Ey/F ux 2Eq 


giving 


Ey \? Eyx Eyx\ |* 
Nu=éuNyx when ¢y= {1-( x) E ais: In ( zy 


Taking all processes into account we have the minimum value given by 
Ni=N yx—(1—46) Nux—(1—¢u)N ux 
as long as H,; >LHy. 
The limiting conditions are satisfied by the above relations: H;=E yy, 
Ey=0, NE=Nuy, HE NZ = 0; Bye Nae er a 


The value of Ey. 

We are assuming that neutral particles do not lose energy by electronic 
excitation. It is known that electronic energy transfer can occur between 
neutral particles, for example the loss of energy of excited Hg atoms in 
the presence of Na atoms, but the cross section is small except when the 
energy transfer is small (Mott and Massey 1949). 

The rate of energy loss by electronic excitation for particles moving 
with velocity v greater than the velocities wu of the valence electrons in 
the solid, is much greater than the rate due to elastic collisions for the 
corresponding velocities, Kinchin and Pease (1955). We shall assume 
that as long as the particle is charged all the energy is lost by inelastic 
processes. However, as the energy decreases the net charge decreases 
and, if H; is the energy at which an atom moving through matter can 
be considered neutral, we assume that all the energy in excess of EH, 
is lost by electronic excitation. 

The work of Stephens and Walker (1955) on fast (6-20 mev) N, O and 
Fl ions has shown that the average charge for a given velocity is in 
reasonable agreement with the criterion proposed by Lamb (1940). 
According to this the moving particle of mass M exchanges electrons 
with the atoms of the solid such that the particle loses those electrons 
whose ionization potential is less than 6(=mEH#/M). If J, is the first 
ionization potential we regard atoms for which H<H,' as, on the average, 
neutral (#;’=MIJ,/m). Although in doing this we shall have over- 
estimated the number of charged particles for H>H,', we shall also have 
underestimated the number of charged particles for H<E,’. 

Bohr (1948) has also considered the value of the average charge on 
an atom which is moving through matter. He gives the average charge 
q as: g=(Z'/ER)m/M where Hp is the Rydberg energy and Z the 


——_- , 
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atomic number of the moving particles. We define a second value of 
E,(H;") below which the average charge is less than +1. The values 
of H,” are within a factor of 2 of those of £,’ as shown in table 4, and as 
also shown in the table this uncertainty does not critically affect the 
values of NJ. 


Table 4. The Results of the Damage Calculations 


Primary Ey ev Ey Ey KH," Nux cm-3 
Si 2-6 10° | 2:9 104 4-2*10° | 3:0x105 9:, D 
O 44105 | 11x104 | 40105 | 2:0«105 36 D 
45 D 
ou dy $) " Nf Nf," 
Si 0-51 1-0 1:0 4-,D 4-,D 
O 0-81 0-99 0-91 29 D 26 D 
34 D 31D 


Both H,' and EH,” are greater than the limit suggested by Seitz (1949) 
below which a charged particle ceased to lose energy by electronic 
excitation—Z ,/""—H,/x where H, for an insulator is the forbidden energy 
gap and x was taken arbitrarily to be 8. 


The value of Ky. 

In his calculations Brinkman (1954) took Hy, to be the energy at which 
the mean distance between non-displacement collisions was equal to the 
interatomic spacing. This value was not put forward as an upper limit 
which, owing to the uncertainties of the details of the process, is the 
value required in the present calculations. 

The criteria for the formation of M regions are that there should be 
non-displacement collisions and that they should occur sufficiently 
frequently so that the energy cannot be dissipated by the lattice. We 
have taken Hy, as the energy below which the elastic collisions can be 
described by the rigid atomic sphere approximation (Ly of Kinchin and 
Pease 1955). This is a convenient upper limit because Snyder and 
Neufeld (1954) have shown that for H>L, the predominant elastic 
collision energy loss is in displacement collisions. Below LL, non- 
displacement collisions occur with increasing probability. The mean 
distance between collisions at Ly is J=N a inds (2 ee le 
giving an upper limit /=30 A. The results of the calculations are shown 


in table 4. 
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The results in the table are not seriously affected by the uncertainties 
in H, and a value of 40 ev has been used for the displacement of either 
atom in quartz. The extreme values are, for all displacements in quartz, — 
31D to 45D. Of the total concentration on the average 2/3 will be 
oxygen displacements, giving 


displaced oxygen atoms: 30 to 21D 
displaced silicon atoms: 15 to LOD. 


The results apply to the low dose region where the rate of replacement 
collisions is negligible. The most serious error is probably the specification 
of the neutron dose which may lead to an uncertainty of 50%, in the above 
values. For the given spectrum, however, the maximum and minimum 
estimates of defect concentrations are unexpectedly close. 


4.4. Remarks on the Nature of the C and E Bands 


Both the C and EH bands increase with neutron dose, indicating that 
these bands are associated with atomic defects since an increase is not 
observed with x-irradiation. The C band can be optically bleached so 
that the observed transition must be due to a trapped electron or hole, 
and not to electrons normally in the neighbourhood of the defect. It is 
assumed that the trapped charge originates from the # defect, which then 
gives rise to the H band. The optical bleaching of the C band, which 
occurs even in synthetic crystals, is attributed to the release of charge 
trapped at C and the recombination of this charge with that trapped at FH. 

The initial rapid and preferential optical and thermal bleaching of the C 
band was found only in neutron irradiated crystals. This bleaching 
probably occurs when the recombination can take place without inter- 
mediate retrapping, such as would be expected when C and £ centres 
are close enough together. It is also required that the mean separation 
for which this process is possible should lead to a modification of the 
energy levels of the C centre so as to give an absorption maximum at 
5-4ev. If this explanation is correct a larger effect would be expected, 
and possibly a 5-4ev peak would be observed, in specimens irradiated 
with electrons of energy near to the threshold for atomic displacements, 
thus favouring the production of close pairs. 

In specimens in which A bands are present, illuminating with C light 
leads to bleaching of both C and A bands. Charge is therefore being 
released from the C centres, being recaptured by the A centres thereby 
neutralizing both. Since, following the O’Brien—Pryce model, an electron 
is required to bleach the A centre we assume that the C defect gives rise 
to the C band after trapping an electron. (There is always some 
absorption in the C region arising from the low energy tail of the D band. 
However, the possibility that the D centre is complementary to A is 
ruled out by the fact that increasing the D band wavelength range content 
of the bleaching light does not increase the bleaching rate of A.) 


—— 
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The conclusion is that the C defect gives rise to the C band after trapping 
an electron and that the # defect leads to the # band after trapping a hole. 

Following the results of heat treatment before irradiation, which 
indicated the presence of silicon vacancies in unirradiated crystals, 
we consider that the # defect is not a silicon vacancy since the EZ band 
is not induced in crystalline quartz by ionizing radiation alone. Further, 
on the grounds that the close complementary defects, suggested for the 
explanation of the rapid preferential bleaching, probably correspond to 
either oxygen interstitial-vacancy pairs or silicon interstitial-vacancy 
pairs we reject the possibility that the C defect is an interstitial silicon 
atom. 

Whether a given defect traps an electron or a hole depends to some 
extent on the nature of the binding in quartz. If it were ionic the 
simplest electron traps are O-* vacancies or interstitial Si**. Having 
already rejected the latter this leads to the conclusion that the C centre 
is an electron trapped at an O-* vacancy, and the # centre an electron 
missing from an interstitial O-?. The properties of quartz—for example, 
the high infra-red absorption and the large magneto-optic anomaly, 
suggest that the binding is mainly ionic. Although the net charge on 
each ion may not be +4 or —2, the traps will trap the same sign of charge 
carrier as above, in which case the C and H defects are associated with 
displaced oxygen ions. 

Taking the concentrations of these displacements as 25D we then 
find for the C and £ defects using (6) and (7) —f,w,0-05, fxg, D0-2 
Since x, and 2, cannot be greater than 1 the lower limits to the f values 
are 0-05 and 0-2 respectively. 


4.5. Defects in Unirradiated Crystals 


The relation between the C and A absorption is still not certain.. One 
possibility is that in the unirradiated crystal the charge on two aluminium 
centres is balanced by an oxygen ion vacancy. This would provide an 
explanation of the observation in I (1954) that, for low neutron doses, 
the absorption in both C and A bands varied with Al concentration, 
In addition, since this mechanism gives a concentration of oxygen 
vacancies related to the substitutional Al content and taking into account 
the overlap of the D band in the C band region, a relation of the form 
(10) above may be obtained. Against the explanation, however, is the 
fact that if true it has to be assumed that oxygen vacancies can be 
generated at 950°C in order to account for p, HIi>p, Il in fig: 14—an 
assumption which is unlikely in view of the insulating properties of 
quartz. Further work is required. 

Both the A, and C bands are broad and they do not become narrower 
on cooling to the temperature of liquid air which suggests that in both 
cases the transitions are between localized levels and one of the energy 
bands of the solid. Concerning the A, band it is probable that only the 
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non-bonding orbitals of the oxygen ions of the aluminium defect are in 
the forbidden energy gap. If the valence band corresponds to electrons 
spending most of their time on oxygen ions the upper valence band will 
correspond to non-bonding electrons having orbitals of the type p, 
and sp,. There will be an allowed transition from this band to the 
p, level of the aluminium centre thereby giving a wide absorption band 
(A,) having the s—p,, anisotropy found by Mitchell and Paige (1955). 

The small optical bleaching in the A, band (table 1) suggests that the 
rate at which holes can move through the upper valence band is low. 
On the other hand the rapid bleaching associated with absorption in the 
C band indicates, independently of the conclusion in § 4.1, that this 
absorption corresponds to excitation of an electron from the centre—the 
excitation being to the conduction band thereby giving the wide and 
temperature independent absorption which is observed. 


§ 5. CONCLUSIONS 


Two of the types of point defect produced by the neutron irradiation 
of quartz have been found to give optical absorption bands in the spectral 
range 1-2 to 8-5eyv. One of the bands (C—5-7 ev) has been shown to be 
associated with an electron trapped at the defect, while indirect experi- 
ments suggest that the other (H—7-6 ev) arises when a positive hole is 
trapped at the second defect. The bands occur at about the same 
energies in fused quartz indicating that the energy of the transition is 
determined by the short range atomic arrangement. 

The existence of more extensive regions of damage following neutron 
irradiation is suggested by the shift of the absorption edge in crystalline 
quartz to the absorption limit which is found in unirradiated fused 
quartz. The latter limit is not changed by irradiation. In these regions 
the structure is considered to be amorphous, the phase change occurring 
as a result of local melting. 

There is some evidence that the defects responsible for the C and # 
bands are oxygen vacancies and interstitial oxygen ions respectively. 

An attempt has been made to put upper and lower limits on the 
concentration of point defects by a collision theory calculation and 
from the results the lower limits to the f values have been found to be 
0-05 (C) and 0-2 (£). 

It is found that all the neutron induced damage may be annealed by 
heating at 950°c. 
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§ 1. INTRODUCTION 


Tux radiative capture of protons of less than 1 Mev in beryllium has been 
studied by a number of authors (see Ajzenberg and Lauritsen 1959, 
Lénsjé et al. 1955, Carlson and Nelson 1955) ; there is some disagreement 
on the number of resonances involved. This paper describes some in- 
vestigations of the gamma-rays produced by the radiative capture of | 
protons in beryllium, which show that their yield for protons of less than 
1 Mev is due to two broad resonances ; any other resonances are weak. 
One of these resonances, at 995 kev, is well known (Fowler et al. 1948, 
Hornyak and Coor 1953); the gamma-radiation from the other, at 
330 kev, has been investigated in detail in the present work. 


§ 2. ExcrraTion FUNCTIONS 


The excitation function of the gamma-rays produced has been investi- 
gated at different times in the course of this work. The detector has 
usually been a 14 in. long x 14 in. diameter sodium iodide crystal attached 
to an EMI photomultiplier type 6262. The target, of beryllium metal 
evaporated onto a copper backing, was bombarded with protons from the 
Cavendish Laboratory two million volt H.T. set. With the detector bias 
set so that pulses corresponding to electron energies of greater than 3 Mev 
in the crystal were counted, the only resonances that were seen were the 
well known broad resonance at 995 kev, the well known narrow resonance 
at 1086 kev and a resonance at 330 kev with a width of approximately 
160 kev, corresponding to a state of 1°B at 6-88 Mev, together with weak 
resonances of variable yield due to the !F(p, «, y) reaction, produced in 
fluorine contamination in the beryllium and on beam-collimating stops. 
A typical excitation function over the lower resonance is shown in fig. 1, 
taken with a target of thickness 10 kev (for protons of 1 Mev). (The high 
points at 340, 480 kev are due to the !F(p, «, y) reaction.) This shows 
well the resonance at 330 kev, while at higher proton energies the yield is 
rising towards the 995 kev resonance. This and similar excitation 
functions are quite comparable with the results of Tangen (Tangen 1946, 
Lo6nsj6 et al. 1955) but disagree with Hunt (1952). 


+ Communicated by Dr. D. H. Wilkinson, F.R.S. 
{ Present address : California Institute of Technology. 
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§ 3. Gamma-RADIATION FROM THE 330 Kev RESONANCE 


The gamma-radiation from the 330 kev resonance has been examined 
in detail, using a sodium iodide crystal spectrometer described in § 2. 
A typical pulse-height spectrum produced by the high energy gamma-rays 
is shown in fig. 2. Marked on it are the positions of the peaks that would 
be produced by gamma-rays emitted to each of the first few states of 
*B. Six distinct peaks are definitely seen in the spectrum. Peaks A 
and B are the total energy peaks of gamma-rays of 1-7, 2-15 Mev which 
were examined in more detail in the low energy spectra. The simplest 
assignment of peaks C, D, E, F is to assume that they are the two- and 
one-quantum escape peaks due to gamma-ray transitions to states of 
B at 0-72, 1-74Mev. This was checked by observing the pulse-height 


Fig. 1 
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Gamma-ray excitation function in the reaction *Be(p, y)'°B. 


spectra with two different crystals: (i) the 13 in. long x 1; in. diameter 
cylinder already referred to and (ii) a lin. cube (both systems having 
similar resolutions). As was expected, peaks C and E were relatively 
higher with the smaller crystal, confirming the assignment of these peaks. 
This spectrum also gives evidence for two more gamma-rays. The bump 
at G appears consistently and could be due to the two-quantum escape 
peak of a gamma-ray emitted to the 2-15 Mev state of ™B ; this gamma- 
- ray is confirmed by the observation of low energy gamma-rays de-exciting 
this state. The high energy tail can be compared with the high energy 
tail of a 6-14 Mev gamma-ray produced by the 340 kev resonance in the 
reaction F(p, «, y). It is certain that there are more high energy pulses 
present than can be due to a 6-14 mev gamma-ray; they are quite con- 
sistent with a weak transition to the ground state of 1B. _ These. results 
are consistent with those of Carlson and Nelson (1955). 
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A general pulse-height spectrum showing the low energy gamma-rays 
is shown in fig. 3; it contains evidence for gamma-rays of 0-72, 1-02, 1:43, 
1-7, 2:15 Mev. Different parts of this spectrum have been examined in 
detail and all these gamma-rays have been seen consistently, together with 
a weak gamma-ray of 0-42 Mev which is not resolved in fig. 3. Using 
results from the °Be(d, n, y) reaction (Rasmussen e¢ al. 1949, Shafroth and 
Hanna 1954) all these gamma-rays can be fitted in as transitions between 
low states of !°B (except for the 1-7 Mev gamma-rays which are discussed 
in § 4) as shown here : 


Gamma-ray energy (Mev) 0-420 0-731 1-018 1-426 2-176 
--0:004 | +0-008 | 40-007 | +0-021 | +0-020 
from state (Mev) | 2-15 0-72 1-74 2-15 2:15 
Assignment 
to state (Mev) 1-74 0 0:72 0-72 0 
Relative intensity 0-045 1-00 0-76 0-043 0-062 


The energies and relative intensities are obtained by comparison with 
pulse-height spectra from known gamma-rays, interleaved with the 
°Be(p, y) runs, knowing the stopping cross section of sodium iodide. The 
relative intensities are thought to be accurate to +25%,. 

Two experiments were done to check that the transitions to the first 
three excited states of 1°B come from the same resonance (evidence is 
given in § 5 for thinking that the ground state transition is not resonant at 
330 kev but is due to the tail of the 995 kev resonance). First the high 
energy gamma-rays were examined at four different proton energies 
ranging from 225 kev to 320 kev. No difference could be seen in the 
pulse-height spectra, implying that if there was a change in the 
relative intensity of the 5-1 and 6-1 Mev gamma-rays it was distinctly 
less than 10%. Similarly six runs were done at four proton energies 
ranging from 245 kev to 330 kev with the total energy peaks of the 0-42 
and 0:72 Mev gamma-rays displayed on the kicksorter. By comparing 
the areas of the total energy peaks it was shown that any change in their 
relative intensity was distinctly less than 10%. These two measure- 
ments show that there is no large change in the relative intensities of the 
4-7, 5-1 and 6-1 Mev gamma-rays over quite a large range of proton energy, 
implying that all these three transitions come from the same state of 9B. 

The high energy gamma-rays were also examined at proton energies 
at approximately 20-30 kev intervals between 250 and 550 key. It was 
found that each pulse-height spectrum was made up of the addition of 
two pulse-height spectra in varying amounts; (a) the pulse-height 
spectrum shewn in fig. 2, (b) the pulse height spectrum from a nearly pure 
transition to the ground state of 1°B. As the proton energy is increased, 
spectrum (a) decreases in intensity while spectrum (6) increases. We can 
estimate what proportion of the pulses in each spectrum corresponds to 
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spectrum (a) and so correct the observed total excitation function to 
obtain the excitation function of the gamma-rays corresponding to 
spectrum (a). Such a corrected excitation function is shown in fig. 4, 
where the continuous line gives the observed total excitation function 
while the points are the estimates of the excitation function of the gamma- 
rays corresponding to spectrum (qa). 

Thus it seems definite that the transitions to the first three excited 
states of 1B are all due to the same resonance of the °Be(p, y) reaction. 
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Pulse-height spectrum of high energy gamma-rays from the reaction 
*Be(p, y)!°B at a proton bombarding energy of 310 key. 
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§ 4. Tun 1-7 Mev Gamma-Ray 
It remains to fit the 1-74 Mev gamma-ray (which has an intensity of 
0-15 of the intensity of the 0-72 Mev gamma-ray) into the decay scheme. 
It cannot be fitted in as a transition from the 1-74 mev state to the ground 
state of !°Be as this would imply a magnetic octupole transition, which 


Fig. 3 
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Pulse-height spectrum of low energy gamma-rays from the reaction 
°Be(p, y)°B at a proton bombarding energy of 310 kev. The arrows 
A, B, C, D, E, F show the expected positions of photopeaks corre- 
sponding to gamma-rays of 0-41, 0-72, 1-02, 1-43, 1-74 and 2-15 mev 
respectively. 


has definitely not been seen in the reaction *Be(d, n, y) (Rasmussen ef al. 
1949, Shafroth and Hanna 1954). So we consider that it is a transition 
from the capture state to one of the states of 1°B at 5-11, 5-16 Mev 
checking this by measuring its energy for different proton bombarding 
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energies. Seven runs were made, with the total energy peaks of the 1-43 
1-7, 2-15 Mev gamma-rays displayed on the kicksorter, at proton energies 
ranging from 245 kev to 374 key, the runs being interleaved with radio- 
thorium calibration runs. It was obvious that, as the proton energy 
changed, the 1-7 Mev peak moved with respect to the 1-43, 2-15 Mev peaks 
so this was checked quantitatively by calculating the gamma-ray energy 
from each run by comparison with the 1-592 Mev two-quantum escape 
peak of the 2-614 Mev gamma-ray in the decay of 7°°T]. A least squares 
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Excitation function of 4:7, 5-1 and 6-1 Mev gamma-rays in the reaction 


*Be(p, y)°B. The full line is the excitation function for gamma-rays 
of all energy in the same reaction. 


analysis, assuming Hy=a--bE,, gave : 
By=(1-486-+0-054)-+ (0-893-L0-174)H#, Mev 
confirming that this gamma-ray is emitted by the capture state. Then 
assuming Hy—A-+0-9H, Mev, from all measurements on this energy of 
this gamma-ray we find : 
A=1-473-40-010 Mev. 
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Then from the @Q-values of the reactions *Be(p, «)®Li, 1°B(n, «)?Li, 
6Li(d, p)7Li and the binding energy of the deuteron (taken from Li e¢ al. 
1951, Ajzenberg and Lauritsen 1955) we find : 
(®Be-+ p—1°B)=6-587-+0-012 Mev 

so that the 1-7 Mev gamma-ray corresponds to a transition to a state of 
0B at 5-114+ 0-016 Mev which agrees very well with the state reported 
at 5-111-+-0-005 Mev and is in definite disagreement with the next state 
at an energy of 5-163--0-005 Mev (Jones and Wilkinson 1954). 

It was also noted that the intensity of the 1-7 Mev gamma-ray relative 
to the intensities of the 1-43, 2-15 Mev gamma-rays did not change with 
proton energy, showing that it came from the same state of 1°B as the 4-7, 
5-1, 6-1 Mev gamma-rays. 

Both from the high energy and low energy spectra it was possible to 
obtain the relative intensities of the gamma-rays emitted by the 6-88 Mev 
state : 


Transition to state of !°B at (Mev) | 0-7: 


Relative intensity 40 | 100 | 20 | 25 


§ 5. GROUND State TRANSITION 

As mentioned above the ground state transition observed at 310 kev 
is believed to be due to the tail of the 995 kev resonance, which is well 
known to emit gamma-rays mainly to the ground state of °B, and not 
due to the 330 kev resonance. This is based on three pieces of evidence : 

(i) The excitation functions and the spectra in the proton energy 
region 350-550 kev reported in § 3 which show that the tail of the 993 kev 
resonance merges into the 330 kev resonance. 

(ii) The observations of Lénsjé et al. (1955) that an excitation function, 
with the bias cut out to correspond to 5-6 Mev excitation energy in the 
crystal so that very few pulses due to 4-7, 5:1, 6-1 Mev gamma-rays would 
be recorded, shows very little indeed of the 330 kev resonance. 

(iii) It is possible to calculate the yield of the ground state transition 
that would be expected at 330 kev due to the 995 kev resonance. This 
was done using both a simple resonance denominator and the full single 
level Breit-Wigner formula, taking the penetrabilities for a standard 
interaction radius of 1-45 (A 13+ 4,13) 108 em from tables of Coulomb 
wave functions. These two estimates only differ by a factor of two at 
330 kev and give roughly the intensity of the ground state transition 
observed ; this was found using the excitation functions described in § 6, 
knowing the relative intensity of the ground state transition at 330 kev. 


§ 6. Rapiative WiptTH 
The radiative width was found by a comparison of the yield of the 
330 kev resonance with that of the 995 kev resonance. To do this we used 
a thick walled brass Geiger counter of the type described by Barnes et al. 
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(1952) with a sensitivity proportional to gamma-ray energy. Using this 
as detector excitation functions were run over both resonances so 
obtaining the relative yield of the two resonances. Then, assuming for 
the 995 kev resonance ’y=23 ev (Hornyak and Coor 1953) we get for the 
330 kev resonance : 

assuming small value of proton width: [y=4-8 ev 

assuming large value of proton width : I[y=1-0 ey. 
We favour the value [’y=4:8 ev as this gives better agreement with the 
theoretical value of the radiative width for the transition to the 1-74 Mev 
state. (Then the radiative width is uncertain to a factor of two, princi- 


pally due to uncertainties in the proton width.) This gives for the partial 
radiative widths : 


Transition to state of !°B at (mev)| 0-72 | 1-74 | 2-15 | 5-11 


Ty (ev) 1:0 | 2-6 | 0-52 | 0:65 


§ 7. ANGULAR DISTRIBUTIONS 
The angular distributions of the transitions to the 0-72 and 1:74 Mev 
states were studied. The theoretical angular distributions considered 
are (capture states 0, 2, 3+ are ignored because of the strength of the 
transition to the 0-++ state at 1-74 Mev) : 


Angular momentum | Capture state To state Angular distribution 
of protons (spin, parity) of spin 
s-wave 1,2— el isotropic 
aa, OG 
p-wave 1+ i) 1+1/10x Tag 08 6) 
1 1 Tee ps0 Peon 8) 
—We [20 x 2(COs 
p-wave -+ + 1/20 ates 


a is the proportion of the compound nucleus formed by channel spin=2 
to that formed by channel spin=1. Writing the p-wave angular distri- 
butions 1+4P, cos@ we see that A has opposite signs for the gamma- 
rays to the 0-72, 1-74 Mev states so that there is a possibility of seeing a 
change in the relative intensity of these two gamma-rays with angle. 
The crystal and photomultiplier were mounted with the front face of 
the crystal 4 in. from the beryllium target so that it could be rotated to 
different angles with respect to the proton beam. Another crystal, 2} in. 
long x 13 in. diameter, was mounted with its front face 1 in. from the 
target to monitor the gamma-ray yield. The isotropy of this set up was 
checked by observing the angular distribution of the gamma-rays from 
the 340 kev resonance of the !%F(p, «, y) reaction, a fluorine target re- 
placing the beryllium target. W hen measuring the angular distribution 
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of the B gamma-rays the pulses from the movable detector corre- 
sponding to electron energies in the crystal from 1-3 to 6-8 Mev were 
displayed on a kicksorter. Eight runs were made spread over five angles : 
0°, 30°, 45°, 60°, 90°; in each case the pulse-height spectrum and the 
number of counts were recorded. The spectra showed no obvious vari- 
ation in the relative intensity of the gamma-rays with angle so a more 
careful comparison was made. Each spectrum was split into six regions 
of pulse-height : four being equally spaced about peaks C, D, EH, F in fig. 2, 
the fifth covering the high energy pulses, and the sixth covering the 
trough between peaks B and C. The angular distribution of the number 
of counts in each region divided by the corresponding number of monitor 
counts was then studied and they were all found to be individually iso- 
tropic. From these results it was possible to say that if the angular 
distributions of the 5-1, 6-1 Mev gamma-rays were of the form 1+-a cos?0, 
then for each of them : 
| a |<0-025. 


Also no departure from isotropy was seen for the 1-7, 4-7 Mev gamma-rays 
but no good limit could be put on their departures from isotropy. 

From these results either of two conclusions is possible for the 6-88 Mev 
state of 1°B : 


either (i) s-wave protons forming a state of 1, 2— 
or (ii) p-wave protons forming a state of 1+, 4:2<a<6:3. 


§ 8. CONCLUSIONS 


If we assume that the 330 kev resonance in the °Be(p, y) reaction is the 
same as the 330 kev resonance in the °Be(p, «) and °Be(p, d) reactions 
(Thomas et al. 1949) which is probable as the widths are the same, then 
the first of the two conclusions at the end of §7 is the more probable as the 
second conclusion seems to lead to a value of the reduced proton width 
which is too large. If s-wave protons are assumed the comparatively 
large radiative width for the transition to the 0+ state makes an assign- 
ment of 1— for the 6-88 Mev state necessary. Then the large deuteron 
reduced width makes a 7’'=0 assignment to the 6-88 mev state follow. If 
this is so the transitions to the 1+, 7'=0 states at 0-72, 1-74 Mev should 
be forbidden by the isotopic spin selection rules on electric dipole tran- 
sitions in self conjugate nuclei, contrary to what is observed. This failure 
of these rules implies a 7’=1 impurity in the 6-88 Mev state of approxi- 
mately 20% in intensity, which is distinctly larger than any other such 
impurity in a light nucleus. This is discussed in detail in a further paper 
(Wilkinson and Clegg 1956). 
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SUMMARY 4 
C and C’ are two electromagnetic cavity resonators which communicate 
with each other through an iris. The present note describes a general 
method of calculating the electric and magnetic field configurations and 
the associated eigenvalues of this compound system, in terms of certain 
orthogonal functions of the separate cavities C and C’ formed when the 
iris is closed by a thin partition. 


§ 1. Lyrropuction 


A PROBLEM which frequently arises with cavity resonators is the 
calculation of the effect on the resonant frequency and @ value of the 
operating mode caused by subsidiary cavities which connect with the 
main resonator. Recessesmay be introduced into the cavity walls for 
the purpose of reducing the effects of unwanted modes, or the method of 
introducing or extracting energy from the cavity may require the addition 
of side tubes to accommodate the coupling elements ; such recesses or 
side tubes may be regarded as small cavities which are coupled to the 


I 


main cavity. Another example is the case of a cavity coupled through 
a slot or iris to a waveguide circuit of low loss which effectively acts as 
a second cavity. Some of these cases have been treated by a boundary 
perturbation method (Slater 1950), or by the method of Hahn (Hahn 
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1941). This note describes a method for treating the general problem of 
two cavities coupled by an iris as in the figure. 


§ 2. THEORY 


Consider the separate cavity C formed when the iris I of the figure is 
closed by a thin perfectly conducting partition. The cavity then has 
completely enclosing walls which will be assumed to be perfectly con- 
ducting. It will further be assumed that it is filled with a uniform 
isotropic medium of permittivity « and permeability wu. If the electric 
field is written as e/1/e sin \/(k/eu)t and the magnetic field as 


h//p cos /(k/ew)t, 
where ¢ denotes time, then the electromagnetic field equations and the 
boundary conditions require that 


V xe=vV(k)h, Veet emis tase ne Go! (CL) 

within the cavity and, 
exds=0=h . ds ee es as Oe) 
at the cavity wall S. The solution of these equations gives rise to an 


infinite set of eigenfunctions e,, h;, and eigenvalues k;. e; and h; determine 
the electric and magnetic field configurations, and can be chosen so that 


{ e; . é; dr=3 =| h; . h; dt . : ° . . (3) 
Cc Cc 


where dz is an element of volume and the integrals are taken throughout 
the whole volume C of the cavity. The frequency of oscillation f; of 
the 7’th mode is given by 
/(k;/e 
WS sis (4) 
Let ¢,, h,, and k;’, similarly be eigenfunctions and eigenvalues for 
cavity C’. Let H, H and K similarly refer to the compound cavity 
formed when the partition is removed. The compound cavity is assumed 
to be filled with a uniform isotropic medium of permittivity « and 
permeability ». Everywhere within the compound cavity 
Vx #=V/(4)A, Ve aA ae 8! Wn) oe () 
and at the walls 
Exds=0=H .ds. Pes Pe rn gO) 
Suppose that throughout the volume of cavity CO, sb, satisfies the wave 
equation V%j,=A,; and the condition Vis;.ds=0 on its boundary. 
The ,’s are assumed to form a complete set of scalar functions in terms 
of which an arbitrary function of position throughout C may be expanded. 
Tf it is further assumed that such an expansion represents the function 
on the boundary of the cavity C as well as throughout its volume, then 
it can be shown that the gradient of a scalar function of position through- 
out C@ may be expanded in terms of the gradients of the functions 7;. 
The 7s may be chosen to satisfy orthonormal relations fo bab; dT=8 y;, 
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in which case the functions g,;=1/—1/,/A; . V¢; satisfy similar relations 

cJi-g;47=8,;. The set ; includes a constant term for which the 
corresponding g is to be taken to be zero. It can be shown (Slater 1950) 
that the g,’s are orthogonal to the h,;’s. That is, fgh;.g;d7=0. Let 
ys,’ and g,’ be defined for cavity C’ in a similar way to ys; and g, for cavity C. 
It can be shown (Teichmann 1953) that within the region C of the com- 
pound cavity, H can be expressed in the form 


A= Sa;hytVur 30.2408 eee 


where w is a scalar function of position. Expanding Vw in terms of the 
functions g,;, eqn. (7) becomes 


H= Ya,b- b> Bg; See 


where, as the functions /;, g; are orthonormal, 6 ;=fog;.H dz which 
can be put in the form Jf, H xg; .ds/./K, the final integration being carried 
out over the surface of the cavity C. Since # xds vanishes on the walls 
of the compound cavity the integral need only be taken over the iris I, 
and thus 6 ;= J; Hxg;.ds/\/K. The «,’s can be expressed 


a =+/(K)f, Exh, .ds/K—k, 
(Slater 1950). Equation (8) now becomes 


a ( Exh,.ds Exg;.ds 
= X4(viw], K—k; jor (, VK ae a 


ds is in the direction of the outward normal to cavity C. Asimilarexpression 
for H can be derived in terms of the functions h,’ and g,’, namely 


PY a pee ne Exg,' .ds 
HT oe / ie Re mad UA Dio |) 
Fv, KTk; yace(f, VK yee} ae 


where the negative sign appears because ds is in the direction of the 
inward normal to cavity C’. Equations (9) and (10) are assumed to 
give the correct tangential component of H on I. To find the perturbed 
eigenvalues, it is now necessary to match the tangential components of 
E and H over I. Everywhere on I, 

~ (vs ea) mi r( Vtk) | a) Ms 


t I 


EXg7).d8\ 4, Exg,;.ds 
+ eet readme... an 


Let 7; form a set of two dimensional orthonormal vector functions of 
position over the surface of the iris I in terms of which any vector tangential 
to | can be expanded. Then 

b= D (Sake, 48) 1 hi’ = > (Sah, «1, ds) Ce © 2D) 


8 8 
and 


9= ZW Gi- 48) Me = LS - 9 8) 0, 


By 


- 
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onl. If Jf; Hx», .ds=y, eqn. (11) becomes 
V(K)frh; .n,dsfyh,’ .n ds 
SrES4( Sr ae, Ji “A n 
t z 


ae V(K) Jr hj. ns dsfrh;. nds SY fri - Ns Es fr gi - 1, ds 
Ks, M JK 1 
Jr 9i- 7s 48 Jr 9i +7 As 
ae (Boe ft nf=o. ae ae acd aren Lo) 
If eqn. (13) is multiplied scalarly by each of the y, in turn and integrated 
over the iris I, an infinite set of equations will be obtained. By eliminat- 
ing the unknown y,’s from these equations a secular equation is obtained. 
i K—k; K—k; ng 
Ir9i Ms 48frgs - neds ae Ir9i- 7548 Ir gi - 7448 
VK V/K 
The roots of this equation give the various eigenvalues of the compound 
system. When one of these has been determined, the ratio’s of the 
y,s may be found from the infinite set of equations derived from (13). 
Knowing the y,’s, the «,’s, 8,’s, «;’s and f,’s can be calculated. This 
enables H to be determined. # can also be found from the «;'s and 
a,;’s (Slater 1950). 


a = tee el 4) 
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[Received June 11, 1956] 


EXPERIMENTS have been carried out on the heat conduction in a tube of 
5x10-%em diameter filled with liquid helium II. The tube whose 
closed end was fitted with a thermometer and a heater was thermally 
insulated, and its open end communicated with a bath of liquid helium 
the temperature of which could be kept constant automatically within 
~2x10-° deg. Applying the two-fluid mode] and assuming that the 
superfluid moving towards the heater suffers no dissipation, it can easily 
be shown that for a tube of length / and radius a, the heat resistance should 
be 

5 = (SS) ap SAE fe en 

p 

where 7, is the viscosity of the normal fluid, p the density and S the 
entropy, the bracketed term being the Poiseuille constant. The aim of the 
work was a determination, using eqn. (1), of the temperature and pressure 
dependence of »,,, which will be described elsewhere. The object of the 
present note is to report observations of the onset of friction made with 
this apparatus. 

According to the suggestion of Gorter and Mellink (1949) there is a 
mutual friction between the normal fluid and superfluid with a force per 
unit volume 

F=Ap,p,(v,—0,)" ee Aides SS (2) 
where v, and v,, are the superfluid and normal velocities, A a constant and 
N is supposed to have the value 3. This would lead to a total heat 
resistance of the form 

a 81 1 

on =n (=) p2S2T +BQ'-1 fr lg 9 ek ons (3) 
where B is a function of the temperature and of the dimensions of the 
tube. It is obvious that if (2) were correct for all velocities, (1) would only 
be valid for vanishing values of Q. On the other hand, experiments have 
shown (Bowers and Mendelssohn 1952, Winkel et al. 1955) that, at least 
in narrow channels, no dissipation of this kind occurs until a critical 
velocity v, is reached. Accordingly 


F=Ap,p,[(v,—v,)—v,]%, (v,—v, >.) ve gee a eS (4) 
has been proposed (see Daunt and Smith 1954) as a refinement of (2). 
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A typical result for A7/Q against Q is given in the figure. Taking 
A and N from the data at high velocities, the curve for eqn. (4) is plotted 
as a dotted line. Neither this curve, nor eqn. (3), shown by the broken 
line, represent the data for all values of Q. Instead, two states of thermal 
resistance appear to be possible. For heat inputs up to a critical value, 
no non-linear friction takes place, and eqn. (1) is followed. At the critical 
heat input, the thermal resistance rises at first almost discontinuously, 
and from then onward increases in agreement with eqns. (2) and (3) if 
N is taken as 4. In the terms of the two-fluid model this means that 
dissipation involving the superfluid, possibly by mutual friction with the 
normal component, appears suddenly and at a finite value. This dissipa- 
tion is velocity dependent and extrapolates back to zero for zero flow 
velocity. The actual dependence of the friction on velocity, which in the 
figure and in a number of other observations requires the value 4 for N 
in eqn. (2), becomes more complicated close to the lambda-point. 


ISO 
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Q =} 
DEG-WAT T 


fele) 


E QU.(3), 
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50 
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Thus, instead of the non-linear resistance appearing gradually beyond. 

a certain velocity, the onset of friction in the superfluid resembles the 
case of a superconductor. There, too, the resistance remains zero up to 
the critical current and then changes over discontinuously to the finite 
value characteristic of the normal state. This similarity is further 
underlined by observations on a tube which was slightly obstructed. 
There the onset of friction was found at a higher value of @ than the dis- 
appearance of friction on decreasing the heat current. A number of such 
hysteresis loops in the thermal resistance were observed at different 
temperatures. Instability in the hysteresis loop could be demonstrated by 
tapping the cryostat, in which case the low heat resistance would suddenly 
change over to the high one or vice versa. While no true hysteresis 
occurred in the unobstructed tube, there was & marked instability in the 
transition region which may be of interest in connection with a recent 


suggestion of Feynman (1955). 
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Note added in proof.—Very large hysteresis effects have since been 
observed in apparently unobstructed tubes of largerdiameter. These effects, 
which under some conditions were time-dependent, were also sensitive 
to vibration, and may possibly account for the inconsistencies in the 
reported values of critical velocity. 
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CXV. An Analysis of the Magnetization Processes in Iron Single Crystals 
by an Electrical Method 


By R. Parker 
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ABSTRACT 


Doring’s theory of the anisotropy of electrical resistivity of ferro- 
magnetics is applied to the published magnetoresistance data of single 
erystals of iron. The theory is shown to be in good agreement with 
experiment and the five anisotropy coefficients for iron are evaluated. 
It is suggested that a study of the magnetoresistance as a function of 
magnetization may be used to obtain information concerning the magneti- 
zation processes in ferromagnetic crystals of cubic symmetry. An 
analysis has been made of the measurements reported by Kuwahara on 
single crystal strips of iron. In certain favourable cases the actual 
magnetization process can be inferred from a number of plausible alterna- 
tives by a consideration of the electrical data. For other specimens an 
estimate of the domain vector orientation in the demagnetized state is 
obtained. . The relative advantages of the electrical method are discussed. 


§1. INTRODUCTION 


Tr is well known that the electrical resistivity of a single crystal of a 
ferromagnetic metal depends upon its state of technical magnetization 
in a complex manner. In order to describe this phenomenon quantita- 
tively it is convenient to define magnetoresistance coefficients, Ap/p, 
where Ap represents the change in the resistivity p, which results from 
the magnetization. Déring (1938) has suggested that the magneto- 
resistance coefficients may be explained by assuming that the resistivity 
is a function of the direction cosines of the magnetization and current 
vectors with respect to the crystal axes. On assuming this function to be 
a power series expanded as far as the fourth power, and upon taking into 
consideration the requirements of Ohm’s Law and of cubic crystal sym- 
metry, an equation can be derived for 4p/p in terms of five independent 
anisotropy constants. These constants have been evaluated for the case 
of nickel and the calculated values are in close accord with experimental 
observations. In comparing theory with experiment Doring has confined 
himself to differences in Ap/p of specimens magnetized to saturation along 
different directions. This has the advantage of avoiding assumptions 
concerning the distribution of domains in the demagnetized state. 
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The only other material for which sufficient experimental data are 
available to test Doéring’s hypothesis is 99-5°% pure iron ; measurements 
were made by Webster (1926, 1927) on single crystal rod specimens cut 
along the three principal crystal directions with magnetic fields applied 
along (longitudinal) and perpendicular to (transverse) the direction of 
current flow. In §2 the five magnetoresistance anisotropy constants 
for iron are evaluated and the experimental results found to be in good 
agreement with Déring’s hypothesis. The extension of the latter to 
initially demagnetized specimens is briefly considered in the light of 
recent developments of our knowledge of ferromagnetic domains. 

In subsequent sections it is shown that Doring’s theory may be applied 
to magnetoresistance data below technical saturation ; in some cases an 
analysis of the electrical data can yield information concerning the 
mechanism of the magnetization processes. In certain favourable cases 
the same information cannot be obtained by any other experimental 
method. 

The magnetoresistance of single crystal strips of iron as a function of 
magnetization has been studied in detail by Kuwahara (1954) and 
Kuwahara and Tatsumoto (1955). The conclusions which may be drawn 
from an analysis of the results of these investigators are discussed in § 4. 


§2. THr MAGNETORESISTANCE ANISOTROPY CONSTANTS OF IRON 


When a single crystal of iron is brought from the ideally demagnetized 
state to technical saturation the magnetoresistance coefficient is given by : 


(p/p) sae=hi( do 7B P —$) + hy > % 058 Pst hs » ot Poe 5” 
7 a 


I FU 4, 94% 
+k,(> 048 7+3 SI a 7a 2—4)+ 2k, »3 a 5% ;(1—a P? —a ;7)8 B, 
t 1,5 FU t,j ze r : ‘ (1) 


where «, and f, are the direction cosines of the magnetization and current 
vector, respectively, with reference to the ith cubic crystals axis, and 
where 4,...k; are constants to be determined from experiment. 
Equation (1) differs slightly from the corresponding equation for nickel 
because of the differences in the easy direction of magnetization in the 
two metals in the demagnetized state. 

Equation (1) is now applied in turn to specimens cut along the three 
principal crystal directions for which measurements have been reported 
by Webster. When a demagnetized specimen with a long axis coincident 
with a [100] crystal direction is brought to saturation in a longitudinal 
field the resulting magnetoresistance coefficient is given by 


(Ap/p)i sat = 3(hy +h). se pug si. |= 8 tae, Sil sare (2) 
For the case of transverse magnetization one obtains from eqn. (1) 
(Ap/P)s sot =—4 ki +4 ke—tk,—(3 ks+a5ks) COS 40)0 0 om ta) 


where @, is the angle between the magnetization vector and one of 
the cubic axes perpendicular to the direction of current flow. The 
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corresponding equations for a specimen cut along a [110] crystal direction 
gel (Ap/p)n ot =b Kr t+h ko tthat+th, . . . . . (4) 
and (Ap/p), sat—=41+a, cos 20,+a,cos46, . . . . (5) 


where 0, is the angle between the magnetization vector and the cubic 
axis perpendicular to the direction of current flow and where 


A= —4(8 hy thy—§ kg t 3 kyt¢ kg) 
Ao= —}(ky—ky +3 ky +@ hq) 
a3g—= ie (3 kg+¢4 ky—ks). 
For a specimen cut along a [111] crystal direction eqn. (1) yields 
(Ap/p)u sat=3 Ket ks+§ ks ek eat Ey 
and (Ap/p)s sxx —4} hott Kz. LE, ee A 


Unfortunately, values for (4p/p)y g4,—(4p/p)1 sa for individual speci- 
mens have not been reported by Webster. However, the reproducibility 
of separate saturation magnetoresistance constants on crystallographi- 
cally similar specimens from Webster’s own results as well as those of 
other investigators (Gondo and Funatogawa 1952) is sufficiently close to 
establish the coefficients to 1x 10-+ in most cases. An exception is the 
value of (4p/p)) «a, for specimens with an easy direction of magnetization 
parallel to the specimen axis. This is undoubtedly due to the variation 
of domain structure in the demagnetized state from one specimen to the 
next and will be considered in greater detail at the end of the section. 
In the meantime no use is made in determining /,...; of any data 
which involve assumptions concerning the demagnetized state of these 
specimens. This still leaves seven independent equations for the evalua- 
tion of only five constants, namely three equations from (5) and one each 
from (4), (6), (7) and one equation from the periodic term in (3). 

After the usual linear extrapolation of the experimental results to 
correct for the changes in spontaneous magnetization (see, for example, 
Bates 1946) one obtains from Webster’s results : 

ky=—8X10-*, k.+57X 10, kgs=+14x Ls, 
Ue as elegee nl OX 10 Cs ae le es (8) 

Saturation magnetoresistance coefficients calculated with the above 
constants and with eqns. (4), (6) and (7) are compared with experimental 
results in the table. A similar comparison is made for eqns. (3) and (5) 
in fig. 1 (a) and 1 (6). It will be noted that the agreement is remarkably 
good, both with regard to numerical values as well as functional depen- 
dence upon the angles 6, and 4,7. It may be concluded that Doring’s 
theory of resistivity anisotropy describes the saturation magnetoresistance 
coefficients of iron crystals so far discussed to well within the uncertain- 
ties of extrapolation and of probable errors of measurement. 


a as ae 8 a ee ae, 

+ Constants differing substantially from these in eqn. (8) are quoted by 
Bozorth (1951); these constants are not in good agreement with some of the 
experimental results described above. 


4G2 


1136 RB. Parker on an Analysis of the Magnetization Processes 


In the determination of the anisotropy constants it has been assumed 
implicitly that in the demagnetized state the volume of domains are 
distributed equally over the [100], [010] and [001] directions. This 


Fig. 1 (a) 


Saturation magnetoresistance. coefficient of [100] specimen. Curve calculated 
from eqns. (3) and (8) ; points from experiment. 


Fig. 1 (b) 


Saturation magnetoresistance of [110] specimen. Curve calculated from 
eqns. (5) and (8) ; points from experiment. 


Comparison of saturation magnetoresistance coefficients calculated from 
eqns. (4) to (8) with experiment 


whe A ’ . Rai 1 
Specimen Field Pe hsctose aee ecpernenn 


results x 10-4 


[110] fenpitacinn ie #. 32 ane 30 t 34 
[111] longitudinal 39 39 to 41 
[111] transverse —15:5 —14to —17 
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assumption is certainly reasonable for well annealed cylindrical specimens 
with their long axes along [110] and [111] directions. For the case of a 
strain and inclusion free specimen cut along a [100] direction however, 
one would expect that in the absence of an external magnetic field the 
magnetostatic free energy is reduced by an almost complete alignment of ' 
domain vectors parallel and antiparallel to the specimen axis. If a 
field is now applied, magnetization proceeds by 180° Bloch wall movements 
and no change in resistivity takes place. If, however, an electric current 
is flowing, a transverse magnetic field is set up, whose value at the 
surface of the specimens investigated by Webster, is several oersteds. 
There is now a tendency for domains to be aligned parallel and anti- 
parallel to the [010] and [001] axes. It is, therefore, not possible to 
predict the domain distribution (and hence the saturation magneto- 
resistance) in the demagnetized state without detailed information with 
regard to coercivity, current density and the ratio of length to diameter 
of the specimens. 

In general, those domains which are initially aligned perpendicularly 
to the specimen axis contribute a term to the saturation magnetoresistance 
coefficient. This according to eqn. (1) is given by 

(Ap/p)usat=hit § hgth hy . Ry aie Set) 


Thus from (8) and (9) the longitudinal saturation magnetoresistance is 
given by the small term 


(A p/p) isa =O X 5°3 X 10-4 5 ‘ <j % c F3 . e (10) 


where c is the fraction of the volume of [010] and [001] type domains 
between the potential leads of the specimen in the demagnetized state. 
On comparison with experiment the results for all specimens can be 
accounted for by values of c between 0 and 1. 

It is significant that c estimated from the results of Gondo and Funato- 
gawa (1952) and Shirakawa (1940) is considerably less than 0-1, whereas 
that from Webster’s results range from 0-3 to 0-9. The former workers 
used specimens of length to diameter ratio in excess of 50, whereas that 
in Webster’s case is about 15. The electrical measurements thus support 
the view that magnetic domains have a preference for those easy direc- 
tions of magnetization which correspond to the smallest demagnetizing 
factor of the specimen. 

The magnetoresistance of single crystal strips of iron as a function of 
magnetization has been studied recently by Kuwahara (1954) and 
Kuwahara and Tatsumoto (1955). The specimens used had approximate 
dimensions 150 mmx2mmxX0:5mm. ‘he remainder of this paper is 
devoted to discussion of the conclusions which may be drawn from our 
analysis of the results of these investigators. 


§ 3. MAGNETORESISTANCE AND THE IMPURITY OF SPECIMENS 


The specimens used by Kuwahara and Tatsumoto were made from a 
block which contained a certain amount of impurity. The effects of 
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impurity upon the magnetoresistance of polycrystalline iron has been 
discussed in some detail by the author (Parker 1952). It was shown that 
the saturation magnetoresistance coefficient of a dilute solid solution in 
iron differs from that of the pure metal by an amount which is principally 
determined by the increase in the specific resistivity that results from the 
presence of the impurity atoms. The impurities listed by Webster have 
a negligible influence upon the resistivity of iron at room temperature. 
It may therefore be concluded that the magnetoresistance coefficients 
reported: by Webster differ from those of ideally pure iron by only a very 
small amount. In the case of the alloy used by Kuwahara however, the 
relevant impurities are 0-66% aluminium and 0-46% silicon by weight. 
The effect of these elements upon (4p/p)y 44 Of polyerystalline iron is 
predicted explicitly by eqns. (2) and (3) in the paper by Parker (1952). 
The predicted values are consistent with the measured values of Kuwahara 
and Tatsumoto over the temperature range —95°c to +100°c. At room: 
temperature both the theoretical and experimental value of (4p/p)) sat 
is 0-60 times that of pure iron. It will be assumed that the value of each 
of the anisotropy constants of eqn. (2) must be multiplied by the same 
factor. This procedure may be justified by the fact that the value 
of Ap/p for every single crystal specimen at any state of magnetization 
shows exactly the same temperature dependence as the value for the 
polycrystalline alloy. The impurity contribution to 4p/p is strongly 
temperature dependent ; it is concluded therefore, that the reduction of 
the polycrystalline magnetoresistance as a result of impurity is due to a 
proportionate reduction in each constituent crystallite. Consequently 
in the subsequent calculations each anisotropy constant in eqn. (8) is 
_ separately reduced by a factor 0-60. 


$4. ANALYSIS OF THE MAGNETIZATION PROCESS OF SINGLE CRYSTAL 
STRIPS 


Of the results to be discussed in this section there is only one case for 
which a detailed knowledge of magnetization processes is already available. 
For this reason it will be considered first. 


4.1. Specimen plane (100), long axis [011 ] 


The distribution of magnetic domains for this specimen has been 
predicted by Néel (1944) and demonstrated with the aid of the Bitter 
figure technique by Bates (1946), Bates and Neale (1950), Bates and Mee 
(1952) and Williams et al. (1949). In the demagnetized state almost the 
whole specimen volume is occupied by domains with their vectors equally 
distributed along the [001], [00I], [010] and [010] directions. In 
addition very small closure domains at the edge of the specimen are 
directed along the [100] and [100] directions. In a small longitudinal 
magnetic field the volume of [001] and [010] domains diminishes by 
180° boundary wall movements until the bulk of the specimen volume 
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consists of alternating [001] and [010] domains. An increase in the 
applied field then causes a rotation of these domains into the direction of 


the applied field. According to eqn. (1) the saturation magnetoresistance 
coefficient is given by 


(4p/p)it sag (1 —0/ 4) (3 ko+4 k3— 4 ky)+(v/2 5) (A; +ho+t k3+8 ky) 

; (TP) 
where v is the volume initially occupied by closure domains and V the 
total specimen volumes, respectively, between the attached potential 
leads. The ratio v/V is of the order of magnitude d/2L, where L is the 
width of the specimen and d the width of [001] domains. For the speci- 
men under consideration the ratio v/V is so small (Lee 1953) that the 
second term in eqn. (11) contributes a quantity too small to be detected 
with the present sensitivity of electrical measurements. It can therefore 
be ignored in the subsequent discussion. 

The magnetoresistance coefficient below magnetic saturation can also 
be determined. While the magnetization proceeds by 180° domain wall 
movements, (4p/p), is zero. In the higher applied field range each 
domain separately is inclined at angles 4, (90°—@), and 90° with respect 


Fig. 2 


[ico] 


Crystal axes and domain vector of (100) [O11] specimen. 


to the cubic crystal axes (fig. 2). Thus according to eqn. (1) the magneto- 
resistance coefficient of the specimen is given by 
(Ap/p)y= 4 by sin 20-+(4 &y-—zz Ma) sin? 26. (12) 
If / is the intensity of magnetization of the whole specimen and J, its 
saturation value then 
T=I, (sin 0+ cos #)X1V/2 . . - ss (13) 
and therefore from (12) and (13) 
(Ap|p)y= he [(2L2—1,2)/21 2 (ks ha) (ere 2b Pas Ct) 
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Values of (Ap/p), calculated from eqn. (14) are now compared with those 
obtained by experiment in fig. 3. It will be noted that excellent agree- 


ment is obtained. 
Fig. 3 


500 IOOO SOO 
Magnetoresistance of (100) [011] specimen. Curve calculated according to 
eqn. (14). 


Experimental results : © increasing field ; * decreasing field. 


4.2. Specimen plane (110), long axis [110] 

This specimen is of some interest since there are a number of plausible 
alternative modes of magnetization. These will now be described in 
turn and it will be seen that only one mode is consistent with the electrical 
data. 

(a) The domain vectors in the demagnetized state are all directed along 
the only directions of easy magnetization contained in the plane of the 
specimen, i.e. [001] and [001]. If an increasing longitudinal field is 
applied the vectors remain parallel to the (110) plane and rotate towards 
the [110] direction. In this case one obtains separately for each domain 
(fig. 4) : 

0 y= = 1/4/2 COB yy, %g=SiN yy, By=P,=1/4/2, Bg=0, . . (15) 


in Iron Single Crystals by an Electrical Method 1141 


Fig. 4 


Crystal axes and domain vector (110) [110] specimen. 


where y,=cos J/J,. On substituting (15) into (1) one obtains after 
simplification 

(Ap/p)w=% +41 COS yy +a, COS? y, . . . . « (16) 
with 

Uy=Hky the th ky tig katt hs) 

A= Thy thet 3 ks +% Ka) 

Ag= — t6(3 byt hy 2 ks). 


This value of (4p/p), is plotted against J in fig. 5, curve A. 

(b) The demagnetized state is characterized by an equal distribution 
of domain volumes over the [010], [010], [100] and [100] directions. 
A small applied longitudinal field then causes the disappearance of 
[010] and [100] domains by 90° or 180° boundary wall movements. 
There is thus no change in electrical resistivity from crystal symmetry 
considerations until J has reached the value J,/\/2. In larger applied 
fields a rotation of the [010] and [100] domains towards the [110] 
direction takes place in the (001) plane. Let the angle between a domain 
vector and the [110] direction be y,. Then from considerations of crystal 
symmetry one finds separately for each domain that 


%y=cos (45°—yo), %2=COS (Ab taey,, ite OSie 2 oc (47) 
On substituting (17) into (1) one obtains 
(Ap/p)y= tha C08 2ya+(F ky—gq Mg)oOs* 2yy, - - + + (18) 


Since [=I,cos y, a relationship is obtained between (Ap/p), and J, 
and this has been plotted in fig. 5, curve B. 

(c) The domain vectors in the demagnetized state are along the [001] 
and [001] directions as described in case (2) but a small applied field 
causes a redistribution in favour of the two easy directions of magnetiza- 
tion nearest to the specimen axis, namely [100] and [010]. In larger 
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fields the magnetization proceeds exactly as described in (6). The 
magnetoresistance coefficient is thus similar to that derived from eqn. 
(18); there is however an additional term }(k,+4,) arising from the 
90° boundary wall movements which predominate in the initial part of 
the magnetization process. The predicted values of (4p/p), as a function 
of I have been plotted in fig. 5, curve C. It will be seen that this is indeed 
the mode which is consistent with the electrical measurements. 


20 


SOO IO00O IS00 
Magnetoresistance of (110) [110] specimen. For explanation of curve sce 
accompanying text. © Experimental results. 


4.3. Specimen plane (011), long axis [100] 

In this case the long axis coincides with the only easy direction of 
magnetization in the plane of the specimen. The magnetization process 
would then be expected to proceed entirely by 180° Bloch wall movements 
within the specimen volume contained between the potential probes. 
The longitudinal magnetoresistance coefficient must therefore be zero for 
all values of J. This is indeed found to be the case by experiment. 
Avery small departure from zero magnetoresistance for [> 0-95 1 , isreadily 
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explained by a disalignment of the crystal axes with the applied field of 
approximately 2°. 


4.4. Specimen plane (001), long axis [100] 

The directions of easy magnetization of this specimen are both along 
(i.e. [100] and [100 ]) and perpendicular (i.e. [010], [010], [001] and 
[001 |) to the specimen axis. From a consideration of shape and magneto- 
crystalline anisotropy only, the domain vector distribution would be 
entirely along the [100] and [100] directions in the demagnetized state. 
As an applied longitudinal field causes 180° boundary wall movements 
only, there can be no change of resistivity according to eqn. (1). 

On the other hand, the [010] and [010 ] directions also lie in the plane 
of the specimen and are almost as favourable directions in the demag- 
netized state, if the length to width ratio of the specimen is not excessive. 
Moreover, the specimen current sets up a magnetic field perpendicular to 
its direction of flow. This field has an estimated intensity of 0-3 oersted 
at the surface of the specimen under discussion. If this field intensity 
should exceed the transverse demagnetizing field appreciably, the domain 
vectors would be aligned perpendicularly to*the direction of current flow 
in the demagnetized state. In this case (Ap/p)y 44, would be given 
according to eqn. (9) as (4p/p)sas= 1 +3 3 +5 a 

If the field due to the specimen current is comparable with the trans- 
verse demagnetizing field, the possibility of coexistence of longitudinal 
and transverse domains in the demagnetized state must be considered. 
The relative volume occupied by these two types of domains is then. 
expected to be sensitive to small elastic stresses and to the magnetic 
history of the specimen. This appears to be the case in the results reported 
by Kuwahara. The latter are shown in fig. 6 together with the value 
calculated from eqn. (9). In the general case 

(Ap/p)sat= (41 +3 Wess Eada eet Se ot ue) 
where ¢ is the fractional volume of transverse domains in the demagnetized 
state. The experimentally determined values of (4p/p), should then 
follow a curve contained within the area OAB in fig. 6. The slope of 
curve can be shown to be ; 

UApleln ok, +H th R@—e, —- -  » 20) 
where c’ represents the ratio of the reduction of volume of transverse to 
antiparallel domains at any point on the magnetization curve. Equations 
(19) and (20) are not independent since the relationship 


Is 
e=11,{ c dl 
0 


must be satisfied. 
With the aid of eqns. (19) and (20) the reported magnetoresistance 


data allow the following inference of the magnetization processes. After 
cooling from the Curie temperature almost the whole specimen volume 
between the potential probes is occupied by domains whose vectors 
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coincide with the specimen axis. Magnetization in an applied longitu 
field then proceeds by the disappearance of antiparallel domains by 180 
Bloch wall movements. If on the other hand the specimen is 
demagnetized in a diminishing a.c. field, over 90% of the volume is 
occupied by domains with their vectors perpendicular to the specimen 
axis. An increasing longitudinal field in the [100] direction will first 
cause the disappearance of all [100] (antiparallel) domains and then a 
growth of [100] domains at the expense of the ‘ transverse ’ domains 
which had their vectors originally orientated along the [010] and [010] 
axes. 

The small additional increase in magnetoresistance coefficient just below 
technical saturation is undoubtedly due to a slight departure of the 
crystal axes from correct alignment with respect to the applied field 
direction. 

Fig. 6 


= 
Ac), x10 
4 


O SOO IOOCO ISOO 


Magnetoresistance of (001) [100] specimen. Line OA consistent with eqn. (9). 
Points according to experiment. 
© Specimen previously demagnetized by a.c. ; 
x Specimen previously cooled from Curie point. 


4.5. Specimen plane (110), long aais [111 ] 

The directions of easy magnetization for this specimen are all equally 
inclined to the long axis of the specimen at an angle of cos— 1//3. 
For this reason any reorientation of the domain vectors among the 
six easy directions leaves the resistivity of the specimen unchanged. Only 
a rotation of the domain vectors away from the nearest easy directions 
towards the field direction will be detected electrically. The magneto- 
resistance coefficients associated with this process will now be calculated. 

Considerations of crystal symmetry require that each domain separately 
is described by direction cosines 


& =H = [(1—cos")/2]}"?,a,= cosd. . . . . (21) 
Since £,=$,=f3;=1/4/3, one obtains from (1) and (21) after due 
simplification : 
(Ap/p)u=($ bo+s kg)+ v/2/8k, sin 26—(§ k2+$kg) cos 26— 
(3 kgta ks) cos 46+ (44/2)k, sin 2p008 GER Tae y aee rea) 
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The magnetization is given by 


I=I,(cosf+/2sind)//3. 2. 2... (23) 
The values of (4p/p), as a function of I have been calculated from 
eqns. (22) and (23) and are shown in fig. 7, together with the values 
reported from experiment. It will be noted that very good agreement is 
obtained. The only conclusion to be drawn from these results is that the 
reorientation of the domain vectors along the easy directions of magnetiza- 
tion nearest to the applied field is almost complete before rotation away 
from the easy directions commences. 


ford 


/ 


Fig. 


500 1000 1500 


Magnetoresistance of (110) [111] specimen. Curve calculated according to 
eqns. (22) and (23) ; points according to experiment. 


§ 5. CONCLUSION 


A detailed analysis of the magnetoresistance data of Kuwahara and 
Tatsumoto on single crystal strips of iron has been given. It was found 
that in a number of favourable cases the electrical data yield detailed 
information concerning the distribution of magnetic domains ne the 
demagnetized state and the nature of the magnetization processes. Good 
agreement is obtained where this information was already available. 
In other cases the analysis of magnetoresistance data facilitates the choice 
between a number of possible alternatives. With one specimen a study 
of the demagnetized state in relation to the magnetic history was possible, 
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In the past the principal experimental method available for the study of 
domain processes has been the powder pattern technique. The magneto- 
resistance method cannot give the same detailed information concerning 
the size and shape of surface domains as the Bitter pattern technique 
since the former only records the integrated effect over the volume of the 
specimen between the potential probes. On the other hand the electrical 
method has a number of advantages. Results are not masked by un- 
representative surface effects ; measurements are readily extended to 
cover a wide range of temperature, applied stress and specimen shape. 
Thus the two methods provide information about the magnetization 
processes, which is largely complementary. A combination of electrical 
and Bitter figure measurements on the same specimen would therefore 
readily extend our present knowledge of magnetization processes in ferro- 
magnetic single crystals of cubic symmetry. 


ACKNOWLEDGMENTS 


The author is indebted to Professor E. R. Andrew for a number of 
helpful suggestions and to the Council of the Royal Society for a grant 
which has enabled him to undertake this investigation. 


REFERENCES 


Barss, L. F., 1946, Proc. Phys. Soc., 58, 153. 

Bates, L. F., and Mrs, C. D., 1952, Proc. Phys. Soc. A, 65, 129. 

Bartss, L. F., and Neate, F. E., 1950, Proc. Phys. Soc. A, 63, 374. 

Bozortu, R. M., 1951, Ferromagnetism (New York: D. van Nostrand Inc.), 
Vp 68: 

Dorine, W., 1938, Ann. Phys., 32, 259. 

Gonpo, Y., and Funatocawa, Z., 1952, Proc. Phys. Soc., Japan, 7, 41. 

Kuwanara, K., 1954, J. Sci. Hiroshima Univ., 18, 87. 

Kuwanara, K., and Tatsumoto, E., 1955, J. Sct. Hiroshima Univ., 19, 385. 

Leg, E. W., 1953, Proc. Phys. Soc. A, 66, 623. 

Neetu, L., 1944, J. Phys. Radium, 5, 241. 

ParKER, R., 1952, Proc. Phys. Soc. B, 65, 616. 

Surrakawa, Y., 1940, Sci. Rep. Téhoku Imp. Univ., 29, 132. 

Wesster, L. W., 1926, Proc. Roy. Soc. A, 118, 196 ; 1927, Jbid., 114, 611. 

Wiiuiams, H. J., Bozortu, R. M., and SHockLEy, W., 1949, Phys. Rev., 75, 


155. 


pea 


CXVI. Angles Between Partial Magnetizations in Spinel Ferrites 
Containing Manganese Ions of more than One Valency 


By W. P. Osmonp 


Mullard Research Laboratories, Salfords, Surrey + 


[Received June 15, 1956] 


ABSTRACT 


In iron-deficient spinel ferrites containing manganese, but no other 
constituent of variable valency, all the manganese will not be divalent. 
Mn** ions tend to form strong directed semicovalent bonds in octa- 
hedral (6) sites but not in tetrahedral (A) sites. The reverse is true for 
Fe?* ions. The antiparallel A-B indirect exchange interaction should 
thus be weaker when one interacting ion is Mn?+ than when both are Fe®+. 
Any Mn+ ions are likely to be coupled with parallel spins by double 
exchange with Mn?* or Mn‘* ions in the same type of lattice site. The 
magnetic ions in either sublattice can then be considered as divided into 
two groups, one based on Mn?* ions, and the other on Fe*+ ions, with 
different degrees of departure from antiparallelism of their net spins to the 
net magnetization of the other sublattice. The molecular field approxima- 
tion as applied by other authors to single ions is extended to such groups of 
ions, with certain provisos as to the correct interpretation to be given to the 
mathematical symbols used for interaction coefficients and angles between 
spin directions. 


§ 1. INTRODUCTION 


Tue chemical formula for a stoichiometric mixed spinel ferrite may be 
written perfectly generally as X, Y, Fe,O, with prqt+r=3, x and ¥ 
representing metalatoms. For compositions containing the stoichiometric 
amount of iron (represented by 50% (mol) Fe,O,) r=2, p+q=l. With 
only ferric iron present, X and Y must then both be divalent (if we exclude 
ferrites containing monovalent ions suchas Li*). In practice the relative 
iron content is often made less than this, so as to reduce the risk of lowered 
resistivity resulting from the easy electron transfer which is known to 
occur between ferrous and ferric ions in similar lattice sites. In such 
circumstances one of the other two metal constituents X, Y must appear 
in other than the solely divalent state, in order that electrostatic neut- 
rality combined with no oxygen deficiency may be maintained. An 
example of this type of composition is provided by the magnesium 
manganese ferrites Mg ,Mn, Fe,O, frequently used for the small cores 
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exhibiting ‘square ’ hysteresis loops in magnetic storage units. In this 
instance manganese is the element appearing with an average valency in 
excess of 2, and unless Mn3+ or Mn*t ions alone are found in similar lattice 
sites there will be lowered resistivity through electron transfer. However, 
any resultant disadvantage is usually considered to be offset by the im- 
proved sintering found with a higher Mn/Mg ratio than that commonly 
employed in compositions of this type for use in microwave components, 
where the dielectric loss caused by electronic conduction is required to be 
negligible. 

Such ‘square loop’ compositions, with all the loss of oxygen arising 
from the deficiency of Fe,0, made good by enhanced manganese valency, 
can be considered as mixtures of MgFe,O, and (a) MnFe,O, and Mn,0, 
for r/2>p, (b) Mn,O, and MgMn,0, for 1>p>7/2, (c) MgMn,O, and 
Mg,MnO, for 1+q>p>1. The last two compositions are not known 
unequivocally to exist, but the formulae provide convenient symbols for 
representing the required manganese valency (Mn for the former and 
Mn!Y for the latter). Compositions for which p>1+q must contain 
some excess MgO, complete oxygen compensation being impossible (if we 
exclude the possibility of the presence of tetravalent iron). In all this 
reasoning it is assumed that the manufacturing precautions have been 
sufficient to prevent the formation of ferrous iron, all this metal then 
being present as Fe** ions. 

From the known ionic distribution in MnFe,O, (Gorter 1954), and from 
the strong tendency of Mn*+ and Mn** ions respectively to form square 
(dsp?) and octahedral (d?s »*) semicovalent bonds in octahedral co- 
ordination with O?- ions (Goodenough and Loeb 1955, to be referred to 
hereafter as G L), we should expect to find the following types of man- 
ganese ions in the octahedral (B) sites of the spinel lattice in the above 
ranges of composition : (a) Mn?+ and Mn+, (b) Mn+, (c) Mn*+ and Mn**. 
In each of these ranges the remaining magnetic ions in B sites will be Fe**, 
and there will also be some of the latter in the tetrahedral (A) sites where 
they have a moderately strong tendency to form stable tetrahedral (sp*) 
semicovalent bonds (GL). In ranges (a) and (b) there will also be some 
Mn?" ions expected in A sites, again with a tendency to form (s p*) semi- 
covalent bonds, while from the resistivity measurements of Economos 
(1954) it appears likely that in ranges (b) and (c), associated with the 
MgMn,0, constituent, there will be some unstable Mn*+ ions in these sites 
as well, probably behaving like Mn?+ + 2Mn3+ < Mn!*, 

Two theories of indirect exchange have an important bearing on the 
magnetic properties of such ferrites. (i) According to Zener’s (1951) 
theory of double exchange, the Mn** ions in B sites will be ferromagneti- 
cally coupled (parallel spins) with up to an equal number of Mn2+ or Mn4+ 
ions in ranges (a) and (c) respectively. The same is also likely to be true 
in A sites in range (b), if Mn+ ions really occur there, while the possible 
unstable manganese ions in these sites in range (c) should also be ferro- 
magnetically coupled. (ii) According to the theory of semicovalent 
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exchange (GL), the antiparallel (antiferromagnetic) interaction between 
semicovalently bonded Fe*+ (or Mn?+) ionsin A sites and the Fe+ and any 
other predominantly ionically bonded magnetic ions in B sites is much 
stronger than that between the magnetic A ions and semicovalently 
bonded Mn+ or Mn** ions in B sites. 

Combining these two theoretical postulates, we can expect that the spins 
of the Mn** ions in B sites, together with those of any other manganese 
ions ferromagnetically coupled with them by double exchange, are less 
likely to be antiparallel to the spins of their nearest-neighbour Fe*+ ions 
in A sites than are those of the Fe**+ ions in B sites. There should thus 
be a reduction of the net moment of the B sublattice compared with the 
value to be expected if Néel’s (1948) law for simple ferrites were followed 
exactly. This reduction can be conveniently expressed mathematically 
as resulting from an angle between the spin direction of the exchanging 
manganese ions (or of the Mn*' ions alone in (b)) and that of the other 
magnetic ions. In the same way, if manganese ions other than Mn** are 
found in A sites in ranges (b) and (c), the ferromagnetic coupling of double 
exchange may lead to a difference between the spin direction of the 
exchanging ions and that of the other magnetic ions in these sites. The 
consequent reduction in net moment of the A sublattice, relative to that 
expected with all spins parallel, can similarly be expressed as resulting 
from an angle between these two spin directions. 

Tn the molecular field approximation used by Néel (1948) for his basic 
theory of the ferrimagnetism of ferrites, invoking antiparallel inter- 
actions between magnetic ions through the superexchange mechanism, 
the existence of A—A and B—B interactions as well as the more important 
A-B interactions is a necessity. This concept of non-negligible inter- 
actions between ions in the same type of site was extended to explain the 
occurrence of angles in either the A or the B sublattice by Yafet and 
Kittel (1952) for the case of a single kind of magnetic ion, and by Niessen 
(1953, 1954) for the case of two kinds of magnetic ions, one kind appearing 
only in B sites. Although with the semicovalent exchange treatment it 
can be shown (Osmond 1956 a) that such A—A and B—B interactions do not 
occur in the same way, the molecular field approximation remains the only 
convenient mathematical treatment of the problem. In this paper 
Niessen’s method is first applied to grouped ions in B sites, and is later 
extended to similar pairs of groups in both sublattices, with the following 
provisos : 

(i) The A—B interaction coefficients are to be taken as average values 
over a relatively small range of compositions. 

(ii) The A—A and B—B interaction terms are introduced to account for 
different spin directions, without necessarily having their usual physical 
significance. 

(iii) The sublattices into which the main spinel sublattices ae further 
divided are to be considered as convenient mathematical concepts, and 
not necessarily as representing actual inter-penetrating crystal lattices 
on which the ions are ordered. 
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§ 2. ANGLES BETWEEN PARTIAL MAGNETIZATIONS IN 
THE B SUBLATTICE ONLY 


In compositions such as those included in range (a) of the last section 
we expect the spins of the magnetic ions (Fe*+ and Mn?*) in A sites to be 
parallel, as they are in B sites in MnFe,O,, on account of their identical 
electronic configuration (3d5) which here enables them both to form stable 
(sp3) semicovalent bonds. We denote by My, the total saturation 
magnetization of this group of ions in Bohr magnetons (8) per formula 
unit. In the same way we denote by Mz, the total saturation magneti- 
zation of the Fe*+ ions and any non-exchanging Mn? ions in B sites, the 
spins of which will again be expected to be mutually parallel. These two 
groups are therefore of precisely the same type. Finally we denote by 
M,’ the total saturation magnetization of the other magnetic ions in B 
sites, consisting of the Mn+ ions and up to the same number of Mn?* ions, 
all spins being parallel through the ferromagnetic coupling of double 
exchange. There is thus a close analogy between this grouping and the 
simpler case of single ions treated by Niessen (1953) in a general way and 
later (1954) in more detail. 

Following Niessen, but with the proviso (iii) given at the end of §1, 
we now imagine the B sublattice to be further divided into two sublattices 
B, and By, with half of each group of ions located in each sublattice in 
such a way that the two partial magnetizations }M, make equal angles ¢ 
with the direction antiparallel to M, and the two partial magnetizations 
4M,’ make equal angles 4’ with the same direction. This symmetrical 
arrangement is shown in fig. 1. It is clear that the resultant moment of 
the B sublattice as a whole is antiparallel to that of the A sublattice, the 
net moment per formula unit for the ferrite being . 

ng—|M,cos¢+M,’cos¢’—-M,|, . . . . (I) 
which naturally reduces to the simple Néel (1948) form when 6=0=¢’. 

We next need exchange interaction coefficients for the molecular field 
approximation. Again following Niessen, subject to the provisos (i) and 
(ii) given at the end of § 1, we define the following coefficients. 

® proportional to the absolute value of the antiparallel (negative) 
exchange integral between two M ions when one is on an A site 
and one on a B site, 
n’ a similar coefficient for use when one of the two ions is an M’ ion, 
n” a similar coefficent which would apply if both exchanging ions 
were M’ ions. 

y; the absolute value of the ratio of the antiparallel B,—B, or By,—By, 
interaction to the corresponding A—B interaction for all combi- 
nation of groups, 

y2 the similar ratio of the B,—B,, interaction to the A—B interaction. 

A logical extension of the first set of definitions is the assumption nn” =n"? 
or z=nin'=n'/n". Strictly speaking we should have three different 
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coefficients y, and three different y, to take account of the different 
possible combinations of interacting ions. However, only one of the y, 
occurs in molecular field components which involve the angles, so that a 
single coefficient is alone required here. We follow Niessen in justifying 
the use of a single y, by the assumption that the physical differences 
between the different types of ions are sufficiently accounted for by the 
differences between the values n, n’, ”. 

According to the molecular field theory the field h,; giving rise to the 
partial magnetization }Mg in sublattice B; has the direction of this 
magnetization. The contribution to this field of each of the other partial 


Fig. 1 


Sublattic e B 


My Sublattice 


Sublattice B. 


pa 


Symmetrical arrangement of partial magnetizations when two different ion 
groups are found in the B spinel sublattice only. 


magnetizations is antiparallel to the particular partial magnetization 
concerned, and therefore has a component perpendicular to the direction 
of fyy on account of the angles ¢, pe The algebraic sum of these per- 
pendicular components must then vanish. Thus ae 
nM , sin d=y2n3 Mz sin 2h6+y n'$M z’ sin (¢’+¢6)—yin'$M ;' sin (p'—¢). 
Writing z=n/n’ and y=sin ¢'/sin 4 this expression reduces i 
22M 4= 227, Mp cos b+ (yety1) Mz’ cos ¢'+y(v2-v1) Me fee d. (2) 
e components of the molecular field h,;' giving rise to 


Phe th se 
Shae age M,’ in the same sublattice we find the condition 


the partial magnetization } ee 

n'M ,sin ¢’=y1n'3M, sin (f'—) Fy 203M sin (6 +4)+y nr" Siz’ sin 24, 

which, since n’/n”=n/n'=z, leads to 

22M ,=(2/y)(y2—v1) Ma cos $! +2(yot+y1) Mp 608 > 
4H2 


+2y,M,' cos¢’. . (8) 


1152 W. P. Osmond on Angles between 


On account of the symmetry of the arrangement the same two expressions 
are derived by considering the molecular fields giving rise to the two partial 
magnetizations in sublattice Bj;. 

Equating expressions (2) and (3) for 22M 4 we find 


y C08 $(¥9—Y1)(2M p+ yM »')=c08 $'(y2—71)(@MptyMe'), - - (A) 
with possible solutions 

(i) cos d/cos ¢’=1/y=sin ¢/sin ¢’, hence =¢’, 

(iil) yo=Y1- 
Niessen (1954) adopts the first solution, thus inherently preserving the 
physical distinction between the sublattices Br and Br indicated by 


17. However, in the resulting simplification of the identical expres- 
sions (2) and (3) y, disappears and we find 


2M ,=y(zM p+’) cos ¢. os 5 Saas ae 


In the present instance of grouped ions, very possibly with different 
relative reductions of the resolved net moment in the direction anti- 
parallel to the moment of the A sublattice, it seems preferable to adopt 
the second solution d<¢’, y;=y2. This means that the two B sublattices 
are equivalent as regards average mutual interaction effects, which is a 
logical deduction from the basic supposition that the subdivision is a 
convenient mathematical fiction rather than a physical reality. Writing 
y=71=Y>2, each of expressions (2) and (3) reduces to 


2M ,=y(zM, cos ¢+M,' cos¢’), . . . . . (56) 
In the type of ferrite being considered we should expect no angles at all 
if there were no Mn** ions present, i.e. when V,’=0. Iftherefore relation 


(5b) may be considered as a legitimate formula right down to this limit we 
find 

y= ,|/M p=Alp, <. Soe Means een ne 
where A and p are the relative proportions in A and B sites respectively of 
the numbers of equivalent ions (Fe*+ and Mn?*) subject to the strongest 
A-B antiparallel interactions. In other words, the degree of departure 
of the spin directions of either B site group from complete antiparallelism 
to all the A site spins increases as the proportion of the Fe*+ group in B 
sites decreases relative to that in A sites. Although, on account of the 
presence of neutral ions such as Mg?*, the relative variation of the three 
groups will not necessarily be linear, the above statement is equivalent 
to saying that the departure from antiparallelism should increase with 
increasing relative proportion of the Mn** group. The mathematical 
treatment thus predicts a trend which is physically plausible. 

Substituting for y in eqn. (5b) we find 


2M ,(1—cos¢)=Mp cOS@, . +44 ws wt 2 ea) 
whence, from (1), the magnetic moment per formula unit is 
ne=|M pi(1—z) cos o-+z}—My|. 2 . . 2 . os (8) 
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If this net moment is to remain antiparallel to that of the A sublattice 
even when Mz<M,, (1—z) cos 6>1—z, whence z>1, which is the actual 
situation assumed in this discussion. When this is true we find from (8) 


(2—1)(1—cos ¢)=(mg+M,—M,)/Mp,. . . . (94) 
whence, from (7) 
(1—1/z) cos ¢’=(ns+M,—M,)/M,’. . . . . (96) 


Tn all these formulae the values of M,, Mz, My’ and n, are assumed. 
to be those per formula unit at the particular temperature considered. 
The temperature variation of the saturation moment of the different ions 
concerned will not be the same, but at 0°K each group moment will be the 
sum of the products of concentration, g factor and spin factor for its in- 
dividual ions. Conveniently the g factor is 2 for Fe*+ and Mn?* ions 
(Gorter 1954) and approximately so (1.99) for Mn*+ (Jonker and van 
Santen 1950), while the spin factor 4/2 is 5/2 for the two former and 2 for 
the latter. For comparison of the theory with experiment we should 
require accurate measurements of the net magnetic moment extra- 
polated to 0°x, determination of the ionic distribution by x-rays or neutron 
diffraction, and knowledge of the non-divalent manganese content by 
chemical analysis. 


§ 3. ANGLES BETWEEN PARTIAL MAGNETIZATIONS IN BOTH A AND B 
SUBLATTICES 


It was shown in § 1 that in the composition ranges (b) and (c) there is the 
possibility of unstable Mn** ions in A sites, as well as the Fe**, and in (b) 
Mn2*, ions expected there, while in both ranges there will be two groups 
of ions in B sites of the types already considered. In such cases it is 
equally logical to divide the A ions into two groups, one containing the 
Mn+ ions, and the other the Fe*+ ions, and to consider that the two groups 
will be subject to A—B interactions of different strengths. We now have 
two groups of the same two kinds in both A and B sublattices. We can 
thus extend the formal molecular field treatment to both these sublattices, 
symbolically dividing.each into two further sublattices with the sym- 
metrical geometrical arrangement of partial magnetizations shown in fig 2. 
The net magnetic moment per formula unit is directed along the axis of 
symmetry and given by 


ng=|M, cos $+ M;’ cos $'—(M 4 cos ~tM 4’ cos ')|. . (10) 


For the mathematical molecular field expressions we retain the inter- 
action coefficients , »’, n”, and y alone in place of 71, 72, with the same 
provisos as before, and we introduce the coefficient « for the absolute 
value of the ratio of any antiparallel A—A interaction to the corresponding 
A-B interaction, for all combinations of groups. The single ‘ average ’ 
coefficients « and y simplify these expressions in spite of the extra terms 
due to the second group of A ions. The condition that the molecular 
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field h,, giving rise to the partial magnetization }M, in sublattice B, 
should be parallel to this magnetization is now 
(nM, cos s+-n'M ,' cosx') sin 6=n'yM ,’ cos 6’ sin d+-ny3M g sin 2¢, 
which, with n/n’=z and 60, gives 

2M, cos4+M 4' cos f’=y(zM gz cos ¢+ M,’ cos 4’). (11) 
On account of the assumption n’/n”=n/n'’=z we find exactly the same 
quasi-symmetrical relationship by considering the molecular field acting 
on the partial magnetization }M,’ in the same sublattice, and naturally 
also, from symmetry, if we consider either of the partial magnetizations 
in sublattice B,;. Applying now the same treatment to the partial 
magnetizations of the A lattice we arrive by similar steps at the relation 

2M, cos ¢+M,’ cos ¢’=a(zM, cos $+ ,’ cos #’), i MAZ) 
from which in conjunction with (11) we find «=1/y, a result which was 
only to be expected from the somewhat idealized symmetrical mathematical 


Fig. 2 


Symmetrical arrangement of partial magnetizations when two different ion 
groups are found in both A and B spinel sublattices. 


arrangement. Moreover, since in the limit when M ,’=0 we expect =0 
and we have the case studied in § 2, we find for this more general case 

oe] ye IM gov Sissel ake a a 
As in the last section, the physical interpretation to be given to this 
mathematical result is that a relative decrease in the number of ions of 
the Fe** group in either sublattice, or increase in the number of ions of 
the Mn** group in the same sublattice, compared with the number of ions 
of the former type in the other sublattice, leads to fewer of the spins in the 
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first two groups being antiparallel to the net moment of the second sub- 
lattice. This is a trend which it is reasonable to expect over a relatively 
small range of compositions, and it is only for small ranges that we can 
assume the value z of the ratios n/n’ and n’/n” to be virtually constant 
when dealing with grouped ions. 

If we now write r,=M ,'/M, and r,=M,'/M, in eqn. (12) we find 

M,;(z cos +r cos ¢’)=aM 4(z cos +r, cos w’), 

whence, from (13) 

z(cos 6—cos f)=r, cos s’—r, cos’, . . . . (14) 
which is the counterpart of relation (7) in §2. We cannot here combine 
this result with (10) to obtain such moderately straightforward expressions 
as (9a) and (9b) of that section. However, if through any change in 
composition the value of z (cos é—cos 7s) remains virtually unaltered, the 
relative changes in the values of 4’ and ¢’ must depend very much upon 
the initial and final values of 74 and r,. 

In any range of compositions for which these simple theoretical results 
are substantially applicable, relation (13) indicates that for M,>M, the 
reduction in the net A lattice moment may well be greater than that of the 
net B lattice moment, since «>1>y so that % would tend to be larger 
than ¢, and possibly also %’ larger than d’. In these circumstances the 
total net moment might then be greater than the value M,+M,'’— 
(M,+M,') given by the simple Néel formula. The opposite might be 
expected for M,>M,. 

§ 4. CONCLUSION 

In this paper the molecular field treatment of the magnetization of 
simple ferrites by indirect exchange has been applied to more complex 
ferrites having two dissimilar groups of magnetic ions, one of which 
contains Fe?+ ions, in one or both of the spinel sublattices. The particular 
type considered is that in which the second group is based on Mn?* ions, 
for which strong nearest neighbour A—Bsemicovalent exchange interactions 
(Goodenough and Loeb 1955), depending on oxygen p? orbitals directed 
towards tetrahedral A sites, are expected to be virtually non-existent. 
The A—A, B-B interaction coefficients and postulated angles are considered 
as convenient mathematical tools for estimating the results of physical 
differences rather than as having a physical existence themselves. 

There are at present no detailed data available by which to check the 
accuracy of the magnetization trends to be expected from this theoretical 
treatment, and for reasons of space we shall not repeat here the results 
found from the values of maximum flux density at room temperature 
reported by Economos (1954) for a series of magnesium manganese ferrites 
in the composition ranges (a) and ()) of § I. This has been done elsewhere 
(Osmond 1956 b), using the results only of the present theoretical work, 
with satisfactory confirmation of the theoretical trends, considering the 
inaccuracies inherent in the necessary averaging factors and assumed 


ionic distributions. 
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ABSTRACT 

The radio echo meteor technique has been used to investigate the 
variation of wind velocity with altitude in the 80-100 km region. Phase 
changes of 5 deg km~! are found in the 12 hour periodic wind and the 
amplitude of this component also increases by about 1m sec! km“! 
with height. There is also a uniform height gradient of 0-5 m sec"* Tens 
in the prevailing wind. The variations in the periodic wind are considered 
in relation to the theory of atmospheric oscillations. 


! § 1. INTRODUCTION 

Two principal methods are available for the systematic observation of 
winds in the 100 km region above the earth. These are the ionospheric 
echo fading technique employing three spaced receivers, and the radio 
echo observation of drifting meteor trails. Because of the properties 
of ionized meteor trails the second method is especially suited to the 
investigation of the variation of wind velocity with height. The height 
of reflection can be determined to within +1 km, independently of time 
of day, season, or radio wavelength. 

Meteor wind observations have been in progress for a number of years, 
principally at Jodrell Bank (Greenhow 1952, 1954, Greenhow and Neufeld 
1955 a) and Adelaide (Elford and Robertson 1953). Measurements 
have shown the presence of large periodic and prevailing wind components 
in the meteor region. Previous results have been restricted to a rather 
limited height range between 90 and 95 km, although there is evidence 
for quite large wind gradients even over this height difference. Of 
particular interest is the height variation of the solar semi-diurnal wind 
component. Preliminary observations (Greenhow and Neufeld 1955 a) 
have shown that the amplitude and phase of this component may change 
considerably with height, and this paper describes an investigation of the 
variations between 80 and 100 km for the year September 1954—August 
1955. 

§ 2. TECHNIQUE 

The coherent pulse system used for measuring drift velocities of meteor 
trails has been described previously (Greenhow 1954). The heights of 
reflection are estimated from echo duration, using an experimental curve 


+ Communicated by Professor A. C. B. Lovell, F.R.S. 
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relating echo decay time and altitude (Greenhow and Neufeld 1955 b) 
In the calibration experiment meteor heights were measured using a split 
beam technique (Clegg and Davidson 1950). The height of reflection 
can be determined to within +3 km from the rate of decay ; which is 
satisfactory for grouping the observations at 5 km height intervals. 

Five duration groups are used (where echo duration is measured 
from maximum to half maximum amplitude) : approximately 0-015—0-03, 
0:03-0:07, 0:07-0:15, 0-15-0-30 and 0:30-0:60 sec. These correspond to 
height ranges 97-5-102-5, 92:5-97-5, 87-5-92-5, 82:5-87-5 and 78:5—82-5 km. 
Above 102:5km the echo duration is too short to enable an accurate 
measure of the wind velocity to be obtained, although echoes are observed 
up to 110 km. The number of echoes recorded below 80 km is very small, 
and integration of results obtained over a long period will be necessary 
in order to determine the wind velocity below this height. 


§ 3. PRESENTATION OF RESULTS 
(a) Velocity Distributions 

North-South and East-West components of wind velocity are measured 
by directing a beamed aerial alternately in these directions. In the 
first years observations all the individual wind measurements from 
50-500 meteors observed in each hourly interval were grouped together. 
The diurnal variations in wind velocity thus referred to a mean height 
between 90 and 95 km, where most meteors are observed (Greenhow and 
Neufeld 1955 a). 

Between September 1954-August 1955 continuous observations were 
made on 2-3 days each month. The echoes obtained during each hour 
have been divided into the height groups listed in §2. Figure 1 shows an 
example of the velocity distributions obtained in this way. The N-S 
wind component is seen to change from 50 m sec-! towards North at a 
height of 100 km, to 50m sec~4 towards South at 85 km showing the 
ability of the method to resolve differences in wind velocity at 5 km 
height intervals. 

Because of the low echo rate in the afternoon and evening, subdivision 
of the results into groups such as in fig. 1 is not always possible. Thus 
while hour by hour wind variations in the 90 and 95 km height groups 
can usually be made throughout the day, at 85 and 100 km an accurate 
wind determination can often be made only between midnight and 
midday. The echo rate in the 80km height group is not sufficiently 
high to enable the mean wind speed in any hour to be determined from 
a single days observations. 


The diurnal variations in wind velocity will now be considered in 
detail 
(b) Component Plots 
Following previous practice NS and EW components of wind velocity 
are plotted at hourly intervals. Examples of the results obtained for 
four heights are given in figs. 2 and 3. Figure 2 shows the observations 
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for November 18-19, 1955 when it was only possible to measure the wind 
in all four height regions during the night. In the N-S direction 12-hour 
oscillations in the wind are present at all heights. The amplitude of the 
oscillation decreases with decreasing height, falling from approximately 
60 m sec“! at 100 km, to 20 m sec! at 85 km. 


Rigel 


97°-5 —102°5 KM. 


(Ty * ©-025. SEC,) 


92°5 — 97°5 KM. 
5 (T3,3 0-05 SEC.) 


87°5 — 92°5 KM. 


(Ty: 0-10 SEC.) 


5 82°5-87°5 KM. 
(11: 0+20 SEC.) 
2 


+ 80 40 {@) 40 BO = 
=I 
M.SEC 
North-South components of wind velocity 0300-0500 hr. Jan. 12-18, 1955. 


i is showing distributions 
Movements towards North +V¥¢.) Histograms show 
of individual meteor drift velocities at heights between 82-5 and 102-5 km. 


Equally well marked is the height variation in phase of this semi-diurnal 
wind component. Maximum velocity towards North occurs near 
0500 hr at 100 km, 0700 at 95 km, 0900 at 90 km, and 1200 at 85 km. 
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This is a mean phase change of approximately 12° km-1, the oscillation 
leading at greater heights. 

Similar results are obtained in the EW direction, but with a phase lag 
of 90° behind the NS component. Thus maximum East velocity always 
occurs 3 hours later than maximum North velocity, corresponding to 


ee 
+404°% ’ / \ 
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Variation of wind velocity with height. 
(a) North-South components of wind velocity Nov. 18-19, 1954. (Movements 
towards North +-ve.) 


(6) East-West components of wind velocity. Nov. 18-19, 1954. (Movements 
towards East +e.) 


a complete rotation of the wind vector in 12 hours. The mean wind shear 
is again 3 m sec? km}, and the change in phase approximately 12° km-1, 

Similar phase and amplitude variations are observed in the complete 
September day shown in fig. 3, although the changes are smaller in 
magnitude. The September observations are discussed in § 4, when the 
method of analysis is described, 


§ 4. Harmonic ANALYSIS OF DAILY OBSERVATIONS 
(a) Method of Analysis 


Twenty-four hour component plots of the type illustrated in fig. 3 
are analysed to give the first three terms of the Fourier series, 
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Vyg=A,4+B, sin 7/12(7-+4,')+C, sin z/6 (+1,") ee 1) 
Vew=4.+ B, sin 7/12 (7+1,')+C, sin 7/6 (7'+4,"). ret) 
In the NS direction, for example, A, is the mean wind speed during 
24 hours, and gives the prevailing wind component. B, is the amplitude 
of the 24 hr periodic wind component, and ¢,’ its phase. Similarly ©, 
and t,” are the amplitude and phase of the 12 hr periodic component. 


Fig. 3 


100 KM. 


~ 
etsy 


Ae ee 


85 KM. 


September 16-17, 1954. Height variation of North-South and East-West 
wind components. The constant and sum of the first two periodic 
terms of the Fourier series are shown as broken curves. 
100km Vyg=+ 1:4+ 6-6 sin 7/12 (t—2-3)+40-0 sin 7/6 (t—0- 
Vaw=+ 2-2+12-5 sin 7/12 (¢—3-4) +26-5 sin 7/6 (t—3- 

95km Vyg=— 2:6+ 2:8 sina/12 (t —3-6) +-23-0 sin 7/6 (¢—1- 
Vew=+10-7+10-9 sin 7/12 (¢ —4-9)+ 23-0 sin 7/6 (¢—3- 

90km Vyg=+ 0-2+ 5-1 sina/12 (t—0-4)+18-1 sin 7/6 (¢—I- 
Vew=+ 89+ 1:9 sin 7/12 (t—3-0)-++16-9 sin 7/6 (t—4- 

85km Vygs=-+ 18+ 47 sin 7/12 (t —3-8)+15'8 sin 7/6 (¢—2-2 
Vew=+13-8+ 1-6 sin 7/12 (t —3-8)+12-0 sin 7/6 (¢—5-0) 


Pe We wo 
wonwnreor 


) 
) 
) 
) 
) 
) 
) 


T’ is the time in hours measured from 0000 hr. Phase is defined such 
that maximum velocity towards North (or East) occurs when Pt t=06 hr 
for the Ist harmonic and 03 hr for the 2nd harmonic. When observations 
are available for only half a day, the diurnal term is omitted. 
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The results for a complete day in September 1954 are given in fig. 3. 
The smooth curves represent the first three terms of the Fourier series 
Vxs» Vew (eqns. (1), (2)). There is good agreement between the 
experimental and theoretical points, showing that harmonics higher than 
the second are negligible. The amplitudes and phases of the NS and 
EW semi-diurnal wind components shown in fig. 3 are represented on 
harmonic dials in figs. 4 (a) (b). The individual points show the time of 
maximum velocity towards North (or East). The changes in amplitude 
and phase with height are very clear, the variations in both directions 
being similar except for the constant phase difference of 3 hours. 

The 24-hour harmonics are shown on a single harmonic dial in fig. 4 (c). 
There are no regular variations in either amplitude or phase with height, 
and no significant difference between the phases in NS and EW directions, 
A true diurnal rotating vector would, of course, show a 6 hour phase 


Fig. 4 
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12-HOUR HARMONIC 12 HOUR HARMONIC 24-HOUR HARMONIC 
N-S COMPONENT E-W COMPONENT _N’S AND E-W COMPONENTS 
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Harmonic dials showing height variation of semi-diurnal and diurnal wind 
components, Sept. 16-17, 1954. 


difference between maximum North and maximum East velocity. The 
mean diurnal wind has a magnitude of only 5m sec~!, and is probably 
due to irregularities in the prevailing and semi-diurnal wind components. 
For this reason the 24-hour harmonic may be neglected when determining 
the prevailing wind from only 12-hours observations. 
(6) 12-hour Harmonic Dials 

The height variations of the semi-diurnal wind component for the 
year September 1954-August 1955 are summarized by the series of 
harmonic dials in fig. 5. Separate dials are used for each height group 
and for both NS and EW components. Inspection of any individual 
harmonic dial shows that large seasonal variations in the amplitude and 
phase of the 12-hour wind component occur. Smooth lines have been 
drawn through the points for individual months or groups of months. 
An attempt has been made to draw similar curves for each harmonic 
dial, so that the trend of the variations at different heights and in both 
directions may be compared at a glance. Between September and 
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—————E——— | 
NORTH — SOUTH M.SEC™! EAST — WEST 
COMPONENTS COMPONENTS 


North-South and East-West components at heights of 


12-hour harmonic dials. 
85, 90, 95 and 100 km, 
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November when the phase of the 12-hour harmonic undergoes a reversal 
the curve has been discontinued. This variation which has repeated 
itself during 3 successive years is extremely complex, and will be discussed 
in a later paper. North-South components for four height groups are 
shown in the left hand vertical column. Inspection of the smooth curves 
shows that the whole pattern decreases in size at lower heights, because 
of the decrease in amplitude of the wind component. Further, in 
successive harmonic dials the pattern slowly rotates, showing a phase 
change of approximately 90° between the 100 km and 85 km height 
groups. 


Fig. 6 
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(a) Change in phase of semi-diurnal wind component over a height difference of 
10km. -+Y¢ change is phase lead at greater heights. 


(b) Change in amplitude of semi-diurnal wind component over a_height 
difference of 10 km. 
90-100 km: @N-S comp. ; OE-W comp. 
85-95 km: xN-S comp. ; +E-W comp. 

Harmonic dials for EW components are shown in the right hand column, 
and similar variations in amplitude and phase are observed at all heights. 
The EW dials have been rotated through 90°, to take account of the 
3-hour period required for the wind vector to rotate from North to East 
As the NS and EW components of wind velocity determined by the present 
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method give two independent measurements of the 12-hour harmonic, 
the similarities between the two series of harmonic dials gives good 
confirmation of the observed changes in amplitude and phase. 


(c) Variation of Phase of 12-howr Component with Height 

The rate of change of phase with height is found to vary with season. 
In fig. 6 (a) the phase change for a 10 km height difference is plotted for 
each month of the year in order to illustrate this seasonal effect. To make 
full use of the observations the difference in phase between 100-90 km, 
and 95-85 km are measured, for both NS and EW components. This 
method gives four determinations of the phase change over a height range 
of 10 km for each month. In fig. 6 (a) a positive change corresponds to 
a phase lead at greatest heights. A smooth line drawn through the 
points shows that the phase difference per km undergoes a marked seasonal 
change, varying from a mean value of 7° km! in winter to 3° km in 
summer. The average annual value is approximately 5° km-!. 


(d) Variation in Amplitude of the 12-hour Component with Height 

The rate of change in amplitude of the semi-diurnal wind component 
with height also varies with season. In fig. 6 (b) the amplitude change 
for a 10 km height difference is plotted for each month of the year. Once 
again the differences between 100-90 km, and 95-85 km are measured 
for both N-S and E-W components. In winter the magnitude of the 
wind gradient is approximately 1-5msec-'km™, falling to only 
0-4 msec-! km-1in summer. Thus in summer-time the amplitude of the 
semi-diurnal wind component is almost constant, increasing from only 
10 to 15 msec-! over the observed height range, while in winter the 
amplitude is approximately 10msec~? at 85 km increasing to over 
30 m sec~! at 100 km. 

If the winter gradient of 1-5 msec-+km™? is extrapolated to lower 
heights, it would suggest that the amplitude of the periodic wind falls 
to zero at 75km. However, the mean gradient over the 90-100 km 
height range is twice as great as that over the 85-95 km region, the respec- 
tive values being 2-0 and 1:0 msec"! km~?. This suggests that the wind 
gradient is reduced at lower heights, and that the amplitude of the 
semi-diurnal wind component falls off less rapidly below 85 km. This 
view is supported by the summer observations which show an almost 
constant small amplitude at all heights. 

It is interesting to note the similarities between the smooth curves 
showing the seasonal variations in amplitude and phase of the 12-hour 
harmonic. These suggest that small amplitude changes are associated 
with small changes in phase, while large phase changes are accompanied 
by correspondingly bigger amplitude variations. 


(c) Height Variation of Prevailing Wind 


(i) Hast-West Component . . 
Because of the difficulty of presenting the four variables, time, speed, 
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direction and height on a single diagram, the prev ailing wind vector has 
also been resolv aA into NS and EW components. As an example, 
the variations in amplitude of the EW component throughout the year 
at heights of 85 and 100 km are shown in fig. 7 (a). It is interesting to 
observe that si eeueh the wind is towards East in summer and winter 
and towards West in spring and autumn at 100 km, at 85 km the wind 
remains towards East in autumn. There is thus a negative wind gradient 
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(a) Seasonal variation of East-West component of prevailing wind at heights 
of 85 and 100 km. 


(6) Gradient of E-W component of prevailing wind. +e gradient: 
increase in east or decrease in west velocity with height. 
(c) Composite diagram showing sign of E-W and N-S components of prevailing 
wind between 70 and 110 km. (Wind blows towards direction shown.) 


at this time. The wind gradient is shown in greater detail in fig. 7 (b). 
As we are dealing with cne component of a vector it has been necessary 
to define a positive wind gradient as one in which the EW wind component 
increases towards East, or decreases towards West. The point for each 
month has been determined by measuring the difference in wind velocity 
over the two 10 km height ranges 85-95 km and 90-100 km. Cretan 
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up to 1 m sec~! km~! are observed during most of the vear with negative 
gradients in February-March and September—October. 

Figures 7 (a), (b) have been combined in 7 (c) to show qualitatively 
the measured wind directions at 85, 90, 95 and 100km; together with 
extrapolated values at 110, 70km. These estimated directions have 
been deduced in the following manner. In September, for example, 
the wind blows towards East and the amplitude decreases with height 
by 0-5 msec"! km~!, falling to only 2msec-! at 100km. Thus the 
mean gradient-observed over a height range of 15 km has only to continue 
for another 4 km, and the wind component will reverse sign to West. 
This view is supported by the October observations, when the reversal 
actually takes place within the observed height range. Additional 
evidence for changes in wind direction at levels outside those directly 
observed is given by inspection of the results for individual days. Thus 
although from June to August the mean wind is towards East at all 
heights, on some days the wind becomes West in the lower height groups. 
Taken in conjunction with the observed mean decrease of velocity of 
1 msec—! km—!, this suggests reversal near 75 km. 


Fig. 8 


(a) Seasonal variation of North-South component of prevailing wind at heights 
of 85 and 97 km. 


(b) Gradient of N-S component of prevailing wind. 


(ii) North-South Component 

The seasonal variations of the N-S component of prevailing wind, 
at heights of 87:5 km and 97-5 km are shown in fig. 8 (a). The curves 
were obtained by integrating the two upper and two lower height groups. 
At 87-5 km the NS component is approximately zero except from May 
to August, when large South components are observed. The positive 
wind gradients of 0:3 m sec-! km-! shown in 8 (b) reduce the amplitude 
of these South components during summer, and produce small North 
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components at other times. The direction of the NS component has 
been superimposed on the EW components in fig. 7 (c). Thus in July, 
for example, the composite diagram shows that the overall wind vector 
is towards South-East above 85 km. 

It is important to note that the wind gradients displayed in figs. 7 (6) 
and 8 (6) not only indicate changes in amplitude of the wind vector, 
but can arise from a change in wind direction. Because of the small 
amplitude of the NS component during most of the year, the variations in 
the overall wind gradient are very similar to those of the EW component. 


§ 5. SUMMARY AND DISCUSSION 


De 
(a) Comparison with Theory 
(i) Amplitude of 12-hour Harmonic 


The semi-diurnal wind component is caused by a 12-hour pressure 
oscillation of the earth’s atmosphere, as a result of thermal and tidal 
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LOG, ,(10#p/p,) PHASE (DEGREES) 
Variation of amplitude and phase of semi-diurnal pressure oscillation with height. 
Experimental observations are compared with theoretical curves (after Wilkes). 
(«) Amplitude variation, py is atmospheric pressure, p the amplitude of the 

12 hr oscillation. 
(b) Phase variation. 
a . ¢ Q > RAT een 4 ~ 
@ Briggs and Spencer, © Elford and Robertson (Southern hemisphere 
summer), < Johnson. 


effects of the sun. At ground level the amplitude of the oscillation is 
only 0-001 of atmospheric pressure, and the associated wind system has 
a theoretical amplitude of only 0-4msec™!, This periodic component 
is obscured by much larger irregular tropospheric winds. At greater 
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heights, however, the amplitude of the pressure oscillation p is increased 
many times by resonance (Taylor 1936, Pekeris 1937). Figure 9 shows 
the predicted variation in amplitude and phase of the semi-diurnal 
pressure oscillation with height (Wilkes 1949). pp») is mean atmospheric 
pressure. The ratio p/p) increases from 0-001 at the earth’s surface to 
0-1 at a height of 100 km with a node at 30km. The ratio p/p) deduced 
from the observed magnitude of the 12-hour wind component between 
85 and 100 km is also shown in fig. 9 (a). The magnification of the 
semi-diurnal pressure oscillation at 100 km varies between 20 and 80, 
agreeing within a factor of 2 with theory. The observed increase in 
amplitude of the wind component with height shows reasonable agree- 
ment with the theoretical predictions. The mean gradient of 
1msec-! km~! is very nearly that required by the change in ratio of 
ppp, due to the increase in resonant magnification with height. 


(ii) Phase of 12-hour Harmonic 


The predicted variations in phase of the 12-hour pressure oscillations 
are shown in fig. 9 (6). Maximum pressure at the earth’s surface occurs 
at 1000 and 2200 br local time, and this corresponds to -++-10° on the 
arbitrary scale. The pressure oscillation undergoes a phase reversal 
at the 30 km node, and above this height maximum pressure 1s at 0400 
and 1600 hr. A slow variation in phase of 5° km~+ occurs above 80 km 
the oscillation leading at greater heights. 

The wind system associated with the semi-diurnal pressure wave in 
the northern hemisphere is such that the wind vector is directed towards 
West when the pressure is a maximum. Thus the observed times of 
maximum North or West wind components enable the phase of the 
pressure oscillation to be determined. Experimental curves for the 
winter and summer, and the annual mean, are shown in fig. 9 (b). The 
theoretical curve agrees well with the mean experimental observations, 
maximum pressure occurring near 0400 and 1600 hr at 80km. The 
observed phase change of 5° km~ is also in agreement with theory. 

‘As in the case of the amplitude variation, seasonal variations in phase 
are observed at a given height. The exact form of the theoretical curves 
in fig. 9 depends critically upon the temperature distribution in the 
upper atmosphere. Thus, seasonal changes in temperature 1m the high 
atmosphere could cause variations of the type observed. 


(b) Comparison with Other Observations 

Mean wind gradients of + 2-3 m sec-+ and 3-6 m sec"? throughout the 
meteor region have been found by the authors (Greenhow 1954) and by 
Elford and Roberston (1953). It is not clear whether these measurements 
refer to the periodic components or prevailing wind, although in the Gase 
surements of Greenhow some sort of average value was 
Briggs and Spencer (1 954), using the ionospheric echo 
he semi-diurnal wind vector is directed towards 


of the mea 
measured. 
technique, find that t 
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North at 0300 hr. in the northern hemisphere. Assuming a mean height 
of 110 km for their observations, this result fits in well with the measure- 
ments described in the present paper, and suggests that the phase change 
of 5°km~! continues throughout the E region. The prevailing winds 
measured by these workers are also approximately twice as fast as those . 
observed in the meteor region, showing a continuation of the steady wind 
gradient to greater altitudes. 

Observations of the semi-diurnal wind system have been made between 
altitudes of 15-20km using normal meteorological methods (Johnson 
1955). At this height the amplitude of the wind vector is only 0-3 m sec~’, 
and its phase is that required by a pressure oscillation which is in phase 
with the pressure oscillation at the earth’s surface. These observations 
are also included in fig. 9, and agree with the theoretical curves below 
the node at 30 km. 


(c) Relation to Sq Variations 


The variation in phase of the semi-diurnal wind component with height 
is of great interest in consideration of the earth’s Sq variations, and the 
situation of the dynamo region. The quiet day variations in the earth’s 
magnetic field can be explained by an ionospheric wind system of the 
type discussed if the wind vector is directed towards North at approxi- 
mately 1100 and 2300 hr. Thus, it is evident that on the basis of existing 
theory the dynamo region does not lie between 85-110 km. However, 
evidence is now available that between these heights there is a 
progressive change in phase of the semi-diurnal wind component; time 
of North velocity advancing from 0730 to 0300 hr. Extrapolation to 
a height where maximum velocity to North occurs at 2300 hr. indicates 
that the dynamo region is at a height of approximately 135 km in 
temperate latitudes. This is within the range suggested by Baker and 
Martyn (1952) but considerably above the height of 105 km observed 
during rocket flights near the magnetic equator (Singer ef al. 1952). 
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EDITORIAL NOTE 


The Editors have received a letter from Professor Dorfman of the Lenin- 
gradsky Gidrometeorologischesky Institut pointing out in regard to the paper 
by D. F. Evans published in the Philosophical Magazine (1956, 1, 370), that the 
principles used by Mr. Evans for measuring nuclear magnetism had been already 
suggested by him in papers published earlier. Professor Dorfman gives 
references to a paper (in Russian) : Doklady, Akad. Nauk. S.S.S.R., 1947, 57, 
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